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1 Introduction 



These notes cover one of the topics of the class given in the Les Houches Summer School 
"Mathematical statistical physics" in July 2005. The lectures tried to give a summary 
of the recent mathematical results about the long-time behaviour of dynamics of (mean- 
field) spin-glasses and other disordered media (expanding on the review |Ben02j ). We 
have chosen here to restrict the scope of these notes to the dynamics of trap models 
only, but to cover this topic in somewhat more depth. 

Let us begin by setting the stage of this long review about the trap models by 
going back to one of the motivations behind their introduction by Bouchaud Bou92j, 
i.e. dynamics of spin-glasses, which is indeed a very good field to get an accurate sample 
of possible and generic long-time phenomena (as aging, memory, rejuvenation, failure of 
the fluctuation-dissipation theorem; see [BCKMOS] for a global review). 

This class of problems can be roughly described as follows. Let F (a compact metric 
space) be the state space for spins and z/ be a probability measure on F. Typically, in 
the discrete (or Ising) spins context F = { — 1, 1} and v = ((5i -|-5_i)/2. In the continuous 
or soft spin context F = /, a compact interval of the real line, and u(dx) = Z^^e~^^^''dx, 
where U (x) is the "one-body potential" . For each configuration of the spin system, i.e. for 
each X = (xi, ...,Xn) G F"^ one defines a random Hamiltonian, Hj{ function of 

the configuration x and of an exterior source of randomness J, i.e. a random variable 
defined on another probability space. The Gibbs measure at inverse temperature (3 is 
then defined on the configuration space F" by 



The statics problem amounts to understanding the large n behaviour of these measures 
for various classes of random Hamiltonians ( |Tal03j is a recent and beautiful book on the 
mathematical results pertaining to these equilibrium problems). The dynamics question 
consists of understanding the behaviour of Markovian processes on the configuration 
space F", for which the Gibbs measure is invariant and even reversible, in the limit 
of large systems (large n) and long times, either when the randomness J is fixed (the 
quenched case) or when it is averaged (often called the annealed case in the mathemat- 
ics literature, but not in the physics papers). These dynamics are typically Glauber 
dynamics for the discrete spin setting, or Langevin dynamics for continuous spins. 

Defining precisely what we mean here by large system size n and long time t is a very 
important question, and very different results can be expected for various time scales t{n) 
as functions of the size of the system. The very wide range of time (and energy or space) 
scales present in the dynamics of these random media is the main interest and difficulty 
of these questions (in our view). At one end of the spectrum of time scales one could take 
first the limit when n goes to infinity and then t to infinity. This is the shortest possible 
long-time scale, much too short typically to allow any escape from metastable states since 
the energy barriers the system can cross are not allowed to diverge. This short time scale 
is well understood for dynamics of various spin-glass models (and related models) in the 
physics literature, in particular for the paradigmatic Langevin dynamics of spherical 
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p-spin models of spin-glasses, mainly through the equations derived by Cugliandolo and 
Kurchan |CK93j (see also [CHS93 ). The fact that the results given by this short-time 
limit are seen as correct for models with a continuous replica symmetry breaking (like 
the Sherrington-Kirkpatrick model) is one of the many bewildering predictions made by 
the physicists. In the Les Houches lectures we covered some of the recent mathematical 
results about this short-time scale for Langevin dynamics obtained in collaboration with 
Alice Guionnet and Amir Dembo r |BD(;n61 IBDGOll IB(;971 [GmOTj). For lack of space 
and in order to keep a better focus we will not touch this topic here at all. 

On the contrary we will be interested in the other end of the range of time scales, 
i.e. time scales t^(n) depending on the system size in such a way that they allow for 
the escape from the deep metastable states. This is where the introduction of the 
phenomenological trap models by Bouchaud becomes meaningful. We will now try to 
explain the relevance of these models in this setting, although in necessarily rather 
imprecise terms. 

At low temperature the Gibbs measure /i^ j should be essentially carried by a small 
part of the configuration space, the "deep valleys" of the random landscape, i.e. the 
regions of low energy, (the "lumps" of p-spin models for instance). So that the dynamics 
should spend most of the time in these regions, which thus becomes very sticky "attrac- 
tors" or "traps" . The trap models ignore the details of the dynamics inside these sticky 
regions. They only keep the statistics of the height of the barriers that the system must 
cross before leaving these regions, and therefore, the statistics of the trapping times 
(i.e. the times needed to escape them). Moreover, the trap models keep the structure 
of the possible routes from one of these traps to the others. The dynamics is, therefore, 
reduced to its caricature: it lives only on the graph whose vertices are all the relevant 
attractors and the edges are the pairs of communicating attractors (see Section |21 for 
precise definition). 

It is a non-trivial matter to prove that these phenomenological models could be 
of any relevance for the original problems. In fact, in the first introduction |Bou92j of 
this model a large complete graph was supposed to be a good ansatz for the simplest 
model of a mean- field spin-glass, i.e. Derrida's Random Energy Model |iDer81j . Proving 
rigorously that Bouchaud's ansatz or "phenomenological" model was indeed a very good 
approximation for metastability and aging questions is quite del icate, an d was done only 
recently, initially in |HH(;n3al IBB(;n3bj . Later we realised (|bcS3) that this very 
simple Bouchaud ansatz for the REM was in fact even better, in the sense that the 
range of time scales where it is a reliable approximation is very wide. 

We now believe that this relevance is even much wider ( "universal" , if we dare) in the 
following sense: some of the lessons learnt from Bouchaud's picture in the REM should 
be relevant for very wide classes of mean-field spin-glasses in appropriate time scales. 
The art is in choosing these time scales long enough so that the (usually sparse) deep 
traps can be found and thus some trapping can take place, but short enough so that the 
deepest traps are not yet relevant and thus the equilibrium (which is of course heavily 
model-dependent) is not yet sampled by the dynamics. This belief has propped us into 
understanding more deeply the trap models in the larger possible generality of the graph 
structure and of the time scales involved. We have indeed found a very universal picture 
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valid for all the examples we have studied, except for the very particular one-dimensional 
Bouchaud trap model which belongs to another class, as we will see below. 

We want to emphasise here that the dynamics of the mean-field spin-glasses is far 
from being the only motivation that makes the study of trap models worthwhile, see 
e.g. |BB03j for references of applications to fragile glasses, soft glassy and granular 
materials, and pinning of extended defects. 

Let us now describe what these notes contain in more detail. We start by giving in 
Section |2 the definition of the Bouchaud trap model for a general graph and a general 
"depth" landscape. We then study, in Sectional the very specific one-dimensional case 
(the graph here is Z) in much detail. The results and the methods are different from all 
other situations we study. We rely essentially on the scaling limit introduced by Pontes, 
Isopi and Newman in |FIN02j . This scaling limit is an interesting self-similar singular 
diffusion, which gives quite easily results about aging, subaging and the "environment 
seen from the particle" . We then go, in Section HJ to the d-dimensional case and show 
that the essence of the results is pretty insensitive to the dimension (as long as c? > 2, 
with some important subtleties for the most difficult case, i.e. d = 2). In particular we 
also give a scaling limit, quite different from the one-dimensional case. We show that the 
properly-rescaled "internal clock" of the dynamics converges to an a-stable subordinator 
and that the process itself when properly rescaled converges to a "fractional-kinetics" 
type of dynamics ( |Zas02j ) which is simply the time change of a (i-dimensional Brownian 
Motion by an independent process, the inverse of an a-stable subordinator. This process 
is also a self-similar continuous process, but it is no longer Markovian. In fact, this 
scaling-limit result is in some sense a (non-trivial) triviality result. It says that the 
Bouchaud trap model has the same scaling limit as a Continuous Time Random Walk 
a la Montroll- Weiss |MW65j . The aging results are then seen as a direct consequence of 
the generalised arcsine law for stable subordinators. 

This picture (valid for all ci > 2) is also naturally valid for the infinite-dimensional 
(or mean-field) case, i.e. for large complete graphs which we study in Sectional Thus, 
we see that the critical mean-field dimension is 2 (in fact, we do not really guess what 
could happen for dimensions between 1 and 2, but it could be an interesting project 
to look at these models on say deterministic fractals with spectral dimension between 
1 and 2). For large complete graphs, which is a very easy case, we choose to give a 
new proof instead of following the well-established route using renewal arguments (as 
in |BD95j and in |BBG03bJ). This proof is slightly longer but illustrates, in this simple 
context, the strategy that we follow in other more difficult cases. An advantage of this 
line of proof is worth mentioning: we get aging results in longer time scales than usual. 

We then use the intuition we hope to have given in the easy Section 15.11 to explain 
in Section 15.21 a general (universal?) scheme for aging based on the same arguments. 
We isolate the arguments needed for the proof given for the complete graph to work in 
general. This boils down to six technical conditions under very general circumstances. 
We then show how, under usual circumstances, these conditions can be reduced to basic 
potential-theoretic conditions for the standard random walk and random subsets of the 
graph. This general scheme is shown to be applicable not only to the cases we already 
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know (i.e. the Bouchaud model on if- with d > 2, or the large complete graphs), but 
also to dynamics of the Random Energy Model (in a wide range of time scales, shorter 
than the one given in [BBUoSEI, including some above the critical temperature) and 



also to very long time scales for large boxes in finite dimensions (with periodic boundary 
conditions) . 



2 Definition of the Bouchaud trap model 

We define here a general class of reversible Markov chains on graphs, which were intro- 
duced by Bouchaud |Bou92j in order to give an effective model for trapping phenomena. 
The precise definition of these Markov chains necessitates three ingredients: a graph G, 
a trapping landscape r and an extra real parameter a. 

We start with the graph G, G = (V, ^), with the set of vertices V and the set of 
edges £. We suppose that G is non-oriented and connected; G could be finite or infinite. 

We then introduce the trapping landscape r. For each vertex positive 
real number which is referred to as the depth of the trap at x. r can also be seen as a 
(positive) measure on V, 

r = ^rA. (2.1) 
Finally, we define the continuous time Markov chain X{t) on V by its jump rates 



.,„=/-^"""°'-?- 'f<--f>^^. (2.2) 
1 0, otherwise. 

Hence, the generator of the chain is 

L/(x)= ^^y{fiy)~fi^))- (2-3) 

Here, the linear scaling factor u defines a time unit and is irrelevant for the dynamical 
properties. Its value in these notes varies for different graphs, mostly for technical 
convenience. The parameter a G [0, 1] characterises the "symmetry" or "locality" of the 
dynamics; its role will be explained later. The initial state of the process will be also 
specified later, usually we set X{0) = 0, where is an arbitrary fixed vertex of the 
graph. 

Definition 2.1. Given a graph G = {V,S), a trapping landscape r and a real constant 
a G [0, 1], we define the Bouchaud trap model, BTM(G, r, a) as the Markov chain X{t) 
on V whose dynamics is given by ()2.2|) and ()2.3|) . 

In the original introduction of the model |Bou92j . the trapping landscape r is given 
by a non-normalised Gibbs measure 

T = 5^r.5. = J]e-^^^5., (2.4) 
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where /3 > is the inverse temperature and is seen as the energy at x. The rates Wxy 
can be then expressed using the random variables instead of t^, 

W:;oy ^ T^exp {(5{{1 - a)Ex - aEy)}, if {x,y)eS. (2.5) 

It is easy to check that r is a reversible measure for the Markov chain X{t)] the 
detailed balance condition is easily verified: 

^xW^xy = J^r^r^ = TyWyx. (2.6) 

Let us give some intuition about the dynamics of the chain X. Consider the embed- 
ded discrete-time Markov chain Y{n): 

X{t) = Y{n) for all S{n) <t<S{n + 1), (2.7) 

where S{0) — and S{ri) is the time of the n*^ jump of X{-). 

If a = 0, the process X is particularly simple. Its jumping rates Wxy do not depend 
on the depth of the target vertex y. X waits at the vertex x an exponentially distributed 
time with mean Txii'dx)'^, where dx is the degree of x in G. After this time, it jumps to 
one of the neighbours of x chosen uniformly at random. Hence, the embedded discrete- 
time Markov chain Y{n) is a simple random walk on the graph G and X{t) is its time 
change. The a — dynamics is sometimes referred to as Random Hopping Times (RHT) 
dynamics. 

If a 7^ 0, the jumping rates Wxy depend on the target vertex. The process therefore 
does not jump uniformly to all neighbours of x. Y{n) is no longer a simple random walk 
but a kind of discrete Random Walk in Random Environment. To observe the effects 
of a > it is useful to consider a particular relatively deep trap x with much shallower 
neighbours. In this case, as a increases, the mean waiting time at x decreases. On the 
other hand, if X is located at some of the neighbours of x, then it is attracted by the 
deep trap x since Wyx is relatively large. Hence, as a increases, the process X stays at 
X a shorter time, but, after leaving it, it has larger probability to return there. We will 
see later that these two competing phenomena might exactly cancel in the long-time 
behaviour of well-chosen characteristics of the Markov chain. 

We are not interested here in a natural line of questions which would be to find 
the best conditions under which the trapping mechanism is not crucial, and the BTM 
behaves as a simple random walk. On the contrary, we want to see how the trapping 
landscape can have a strong influence on the long time behaviour. Obviously, this can 
happen only if this trapping landscape is strongly inhomogeneous. 

Strong inhomogeneity can be easily achieved in the class of random landscapes with 
heavy tails, the essential hypothesis being that the expectation of the depth should be 
infinite. One of the assumptions we will use is therefore: 

Assumption 2.2. The depths Tx are positive i.i.d. random variables belonging to the do- 
main of the attraction of the totally asymmetric a-stable law with a e (0, 1). This means 
that there exists a slowly varying function L (i.e., for all s > lim„_>oo L{us) / L{u) — 1), 
such that 

¥[tx>u]=u~''L{u). (2.8) 
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Sometimes, to avoid unnecessary technical difficulties, we use the stronger assump- 
tion: 

Assumption 2.3. The depths Tx are positive i.i.d. random variables satisfying 

lim u°P[r^ >u] = l, aE (0, 1). (2.9) 

u— >oo 

These assumptions are satisfied at low temperature for the standard choice of the sta- 
tistical physics literature: for —E^ being an i.i.d. collection of exponentially distributed 
random variables with mean 1. The depth of the traps then satisfies 

P[r^ >u]= Fle'f^^- >u]= u'^^^, u>l. (2.10) 

Hence, if /? > 1, then Assumption 12.31 is satisfied with a = 1/(3, and the expected value 
of the depth diverges. 



2.1 Examples of trap models 

Specifically, we will consider the following models: 

1. The Bouchaud model on Z. Here G is the integer lattice, Z, with the nearest 
neighbour edges. The trapping landscape will be as in Assumption 12.21 We will 



i pe 

report in Section El on work by |KI Nn2[ IRCn^l ICerOfi 



2. The BTM on Z'^, d > 1, with a random l andscape as in Assumpti on 12.31 The main 



results about this case are contained in |RCMn6l ICern3l iRCOfib 



We will also deal with a generalisation of the former setting, i.e., with a sequence of 
Bouchaud trap models, {BTM(G'„, r^, a) : n G N}. We will then consider different time 
scales depending on n and write Xn{t) for the Markov chain on the level n. In this case 
the law of may depend on n. Then the requirement E[r^] = oo is not necessary, as 
we will see. 



3. The BTM in a large b ox in Z*^ with periodic boundary condition, here Gn = Z'^/ nZ" 



is the torus of size n BC06a 



4. The BTM on a large complete graph. Here we consider a sequence of complete 
graphs with n vertices. This is the model that was originally proposed in |Bou92j . 

5. The Random Energy Model (REM) dynamics. We deal here with a sequence of n- 
dimensional hypercubes Gn = { — 1, The landscape will be given by normally 
distributed (centred, with variance n) energies E^. This is the case where r's are 
not heavy-tailed. 

We will see that all the previous cases with exception of the one-dimensional lattice 
behave very similarly. 
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2.2 Natural questions on trap models 

We will be interested in the long-time behaviour of the BTM. There are several natural 
questions to ask in order to quantify the influence of trapping. The most important 
question we will address in these notes is the question of aging (at different time scales). 
Let us now define what we will call aging and subaging here. We will consider a time 
interval [t^, + 1], where the waiting time tyj (or the age of the system) as well as the 
length t of the time window (or duration of the observation) will grow to infinity. We 
will then consider various two-time functions, say C{tuj,tw + t), which depend on the 
trajectory of the Bouchaud trap model in the time window [tu,, tu, + t]. We will say that 
there is aging for these two-time functions iff 

lim C{t^,t^ + 9t^) (2.11) 

exists and is non-trivial. We will call this limit, C{9), the aging function. We will say 
that there is subaging with exponent 7 < 1 iff 

lim C{t^,t^ + 9tl) (2.12) 

tyj 

exists and is again non-trivial. 

We need now to define good two-time functions in order to be able to deal with aging 
for BTM. The following functions are mostly studied: 

(a) The probability that, conditionally on r, the process does not jump during the 
specified time interval [tw, tw + t], 

U{tw,t^ + t;T) =F[X{t') = X{tw)yt' e [t^,t^ + t]\T]. (2.13) 

(b) The probability that the system is in the same trap at both times t^ and tw + 1, 

Ritw,tw + t;T) ^F[X{tw) ^ X{tw + t)\T]. (2.14) 

(c) And finally, the quantity 

i?«(i^,i^ + t;T) =E[^[P(X(i^ + i) =x|T,X(i^))]2|T], (2.15) 

which is the probability that two independent walkers will be at the same site after 
time t + tw if they were at the same site at time tw, averaged over the distribution 
of the common starting point X{tw)- 

These quantities are random objects: they still depend on the randomness of the 
trapping landscape r. They are usually called guenc/ieo? two-time functions. In addition 
to the quenched two-time functions, we will also consider their average over the random 
landscape. We define the averaged two-time functions: 

n(t^, tw + t)^ F[x{t') = x{tw) yf e [tw, tw + 1]] , 

R(tw, t + tw)^ nx(tw) = x(tw + 1)], ^2 16) 

+ t)^E [j2[¥{X{tw + t)^ x\t, X{tw))]''] . 
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We will state aging results for both quenched and averaged two-time functions for 
the various Bouchaud trap models given in Section ITTl We will strive to get the widest 
possible range of time scales where aging occurs. 

Even though our main motivation was the study of aging, there are many other 
questions about the long-time behaviour of the BTM which are of interest. For instance: 

• The behaviour of the environment seen from the particle. Here the prominent 
feature of this environment seen from the position at time t is simply the depth 
of the trap Tx(t) where the process is located. We will give limit theorems for this 
quantity, which are crucial for most of the aging results. 

• Nature of the spectrum of the Markov chain close to its edge. Naturally, the long 
time behaviour of X(t) can be understood from the edge of the spectrum of the 
generator L. This question deserves further study (see jBF05| IBFOGj and also 
MB97p . We intend to address this question for BTM in finite dimensions in a 
forthcoming work. 

• Anomalous diffusion. In the case where graph is Z*^, can we see that X(t) is slow: 
for instance that E[X(t)2] < t? Or get the tail behaviour of \X(t)\7 

• Scaling limit. Again in the case of Z'^, is there a scaling limit for the process X{t), 
i.e. a way to normalise space and time so that X'^{t) = eX(t/h{e)) converges to a 
process on which we can describe. 

2.3 References 

The physics literature on trap models is so abundant that we cannot try to be exhaus- 
tive. For earlier references on finite-dimensional questions see |Mac85| IBGQOj where the 
anomalous character of the diffusion is given as well as a scaling limit. For aging ques- 
tions on large complete graphs and relation to spin-glass dynamics see fBou92. ,BD95j . see 
also |B(yKM9'8] for a more global picture. For aging questions for the finite-dimensional 
model see |MB961 IHMBOOL IRMBOIL IBBOaj among many other studies. 

We will give references to mathematical papers at the end of every section, when 
necessary. 

3 The one-dimensional trap model 

We will consider in this section BTM(Z, r, a) on the one- dimensional lattice, Z, with 
nearest-neighbour edges. The random depths, r^^, will be taken to be i.i.d., in the 
domain of the attraction of an a-stable law, a < 1, as in Assumption 12.21 There are 
several reasons why this particular graph should be treated apart. First, as usual, the 
one-dimensional model is easier to study. Second, as we have already mentioned, the 
one-dimensional BTM has some specific features that distinguish it from all other cases 
presented later in these notes. 
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Another distinguishing feature (of technical character) is that, at present, the one- 
dimensional BTM is the only case where the asymmetric variant (a > 0) has been 
rigorously studied. For technical convenience we choose here v = Va = 1E[tq"'*]2/2, that 
is we set 

w,y = iE[ro-"]V-(^-"V», if \x-y\ = 1. (3.1) 

We set X{0) = 0. We suppose that z/^ is finite for all a G [0, 1], this is obviously the case 
if, e.g.. To > c a.s. for some c > 0. 



3.1 The Fontes-Isopi-Newman singular diffusion 

The most useful feature of the one- dimensional BTM is that we can identify its scal- 
ing limit as an interesting one- dimensional singular diffusion in random environment 
introduced by Pontes, Isopi and Newman jFINn2> . 

Definition 3.1 (The F.I.N, diffusion). Let {xi,Vi) be an inhomogeneous Poisson 
point process on M x (0, cxd) with intensity measure dx av~^~°'dv. Define the random 
discrete measure p = Yli'^i^xf We call p the random environment. Conditionally on p, 
we define the F.I.N, diffusion Z{s) as a diffusion process (with Z{0) = 0) that can be 
expressed as a time change of a standard one-dimensional Brownian motion B{t) with 
the speed measure p, as follows |IM65j : Denoting by £{t, y) the local time of the standard 
Brownian motion B{t) at y, we define 



Ut)= / ^{t,y)p{dy) (3.2) 

and its generalised right-continuous inverse 

^p{s) = m{{t > : (j)p{t) > s}. (3.3) 

Then Z{s) = B{iljp{s)). 

The following proposition lists some of the properties of the diffusion Z and the 
measure p that may be of interest. 

Proposition 3.2. (i) The intensity measure of the Poisson point process is non- 
integrable at v = 0, therefore the set of all atoms of p is a.s. dense in M. 

(a) Conditionally on p, the distribution of Z{t), t > 0, is a discrete probability measure 
vl = J2i'^i(f)^xi, with the same set of atoms as p. 

(Hi) The diffusion Z has continuous sample paths. 

(iv) Define Pp{t,Xi) = F[Z{t) = Xi\p]/vi for all atoms Xi of p and t > 0. The function 
Pp has a unique jointly continuous extension to (0, oo) x M. Moreover, pp satisfies 
the following equation 

|pp(t,x) = ^p,(t,x). (3.4) 
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The singular differential operator d"^ f / dpdx is defined (see, e.g., \DM7(\ \KW8^ ) 
by h = d'^f/dpdx if for some Ci, C2 G M 

(c2 + J h{v)p{dv)'^du. (3.5) 

(v) The diffusion Z and its speed measure p are self-similar: for all X > 0, t > and 

p([0,x]) A-^/V([0,Ax]) and Z(t) A-^Z(tA(^+°)/"). (3.6) 
Therefore, the diffusion Z is anomalous. 

(vi) There exist constants C , c such that for all x and t > 

¥\\Z{t)\>x] < Cexp [-cf^V'^'T (3-7) 

Proof. Statement (i) is trivial, (ii) is proved in |FIN02j . Claim (iii) follows from (i), 
the continuity of sample paths of B, and the definition of Z. (iv) is a non-trivial claim 
of the theory of quasi-diffusions, see above references. The first part of (v) is a direct 
consequence of the definition of p. The second part then follows from the first one and 
from the well-known scaling relations for the Brownian motion B and its local time: 

B{t)'= X~^B{XH) and i{t,x)'= X-H{XH,Xx). (3.8) 



The last claim is proved in jCerOGj . □ 

Remark. The scale of the upper bound in (vi) is probably optimal. The corresponding 
lower bound was however never proved. The numerical simulations and non-rigorous 
arguments in |BB03j however support this conjecture, and give even exact values for 
constants C and c. 

3.2 The scaling limit 

We now explain how the F.I.N, diffusion appears as a scaling limit of the BTM. For all 
e G (0, 1) we consider the rescaled process 

X'{t)=eX{t/ece), (3.9) 

where 

c, = (inf[t > : P(ro > t) < e])~\ (3.10) 

It follows from Assumption 12.21 that there is a slowly varying function L' such that 
Cg = e^^"L'(l/e). We will also consider the following rescaled landscapes, 

T^(dx) = c,5^r,4.(dx) =: ^r44.(dx). (3.11) 



It is not so difficult to see that the distribution of converges to the distribution 
of p as £ ^ 0. The next proposition states that it is possible to construct a coupling 
between different scales such that the convergence becomes almost sure. 
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Proposition 3.3 (Existence of coupling). There exists a family of measures 
and processes X'^ constructed on the same probability space as the measure p and the 
Brownian motion B such that 

(i) For all e > 0, has the same distribution as r^. 

(ii) The measures converge to p vaguely and in the point-process sense, p-a.s. 

(Hi) can be expressed as a time(-scale) change of the Brownian motion B with the 
speed measure (see Section \S.2.1\) . It has the same distribution as X^. 

We can now state the principal theorem of this section. It was proved in |FlNn2j for 
a = and in [BCOSi for a > 0. 

Theorem 3.4 (Scaling limit of the one-dimensional BTM). As e 0, for every 
fixed t > and all a E [0, 1], the distribution of f^^^^-j) converges weakly and in 

the point-process sense to the distribution of (Z(t), p({Z(t)})) , p-a.s. 

The notion of convergence in the point-process sense used in the theorem was in- 
troduced in |FINn2j . It is used here because we want to deal with quantities like 
P[X^(t) = X^{t')] to prove aging. The usual weak or vague convergences of measures 
are insensitive to such kind of quantities. This notion of convergence is defined by 

Definition 3.5 (Point-process convergence). Given a family u, z/^, e > 0, of locally 

finite measures on M, we say that z/^ converges in the point process sense to z/, and write 
^ as e ^ 0, provided the following holds: If the atoms of z/, u'^ are, respectively, at 
the distinct locations yi, yfi with weights Wj, then the subsets of [/^ = Uj'{(|/^/, wf/)} 
of M X (0, oo) converge to U = Wj)} as e — in the sense that for any open 

O, whose closure is a compact subset of M x (0, oo) such that its boundary contains no 
points of U , the number of points It/*^ fl 0| in f/^ fl O is finite and equals |f/ fl 0| for all 
e small enough. 

Remark. The convergence of ^^x'^^^^ gives a description of the environment seen by the 
particle. More explicitly, the distribution of the normalised depth of the trap where X 
is located at large time CeTx{t/ece) converges to the distribution of p{{Z{t)}). 

We now sketch the three main tools that are used in the proof of Theorem 
3.2.1 Time-scale change of Brownian motion 

To better understand how (Z, p) arises as the scaling limit of (X, r), one should use the 
fact that not only diffusions, but also nearest-neighbour random walks in dimension one, 
can be expressed as time(-scale) change of the Brownian motion. The scale change is 
necessary only if a 7^ 0, because the process X{t) does not jump left or right with equal 
probabilities. 

We first define the time-scale change. Consider a locally-finite, discrete, non-random 
measure 




(3.12) 
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which has atoms with weights Wi at positions y,. The measure will be referred to as 
the speed measure. Let S' be a strictly increasing function defined on the set {?/«}. We 
call such 5* the scaling function. Let us introduce slightly nonstandard notation S o fi 
for the "scaled measure" 

{S o/i)(dx) = ^Wi5s{y,){dx). (3.13) 

i 

Similarly as in definition of Z, we define the function 

0(/i,^)(t)= [ e{t,y){So^^){dy) (3.14) 



and the stopping time ■ip{fi,S){s) as the first time when (f){fj,, S)(t) = s. The function 
0(/i, S){t) is a nondecreasing, continuous function, and il^{fj., S){s) is its generalised right 
continuous inverse. It is an easy corollary of the results of |Sto63j that the process 

X{fi,Sm:=S-\B{^{f,,S)m (3.15) 

is a nearest-neighbour random walk on the set of atoms of fi. Moreover, every nearest- 
neighbour random walk on a countable, nowhere-dense subset of R satisfying some mild 
conditions on transition probabilities can be expressed in this way. We call the process 
X{fi, S) the time-scale change of the Brownian motion. If = Id, the identity mapping, 
we speak only about the time change. 

The following proposition summarises the properties of X{^,S) if the set of atoms 
of fi has no accumulation point. In this case we can suppose that the locations of atoms 
yi satisfy yt < yj if i < j. 

Proposition 3.6 (Stone, [Sto63j). The process X{fi, S){t) is a nearest-neighbour ran- 
dom walk on the set {yi} of atoms of fi. The waiting time in the state yi is exponentially 
distributed with mean 

o,„ iS{yi+i) - S{yi)){S{yi) - Sjyi^i)) 

5(y.^0 - S{y..,) • ^^-'^^ 

After leaving state yi, X(/i, 5*) enters states and with respective probabilities 

S{yi+i) - S{yi) S{yi) - g(y»-i) 

and — — -. (-J-l') 



SiVi+i) - S{y^-l) S{yi+i) - S{yi-i)' 

Using this proposition it is possible to express the processes X'^(t) (see ()3.9j] ) as a 
time-scale change of the Brownian motion B. It is not surprise that should be chosen 
as the speed measures. The scaling function is defined by 

Six) = i - °' (3.18) 

otherwise. 
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where 

r. = l^'r-^r-^,. (3.19) 

Observe that Ua was chosen in such way that ¥\rx] = 1. If a = 0, S' is the identity 
mapping on Z, there is no scale change in this case. Define further S^{-) = eS{e~^-). 
It is easy to check, using Proposition I3.fj| that the processes X{t'^,S'^) have the same 
distribution as X^. 

It is convenient to introduce processes X^{t) that are the time change of the Brownian 
motion with speed measures S"^ o r^. Namely, 

X'{t) =X{S' oT',ld){t). (3.20) 

The processes are related to X" by X%t) = {S^)-\X%t)). 



3.2.2 Convergence of the fixed-time distributions 

We have expressed the processes Af^ as the time change of the Brownian motion with 
the speed measure S*^ o r^. We want to show that X'^ and mainly converge to Z. 
As stated in the following important theorem, it is sufficient to check the convergence of 
the speed measures to prove the convergence of fixed time distributions. Observe that 
the theorem deals only with non-random measures. 

Theorem 3.7 ([Sto63, FIN02j). Let fi'^ , fi be a collection o/ non-random locally-finite 
measures, and let , y be defined by 

y%t) = X{fi', Id) (t) and y{t) =X{fi,ld){t). (3.21) 

For any deterministic to > 0, let z/^ denote the distribution ofy^lto) and v denote the 
distribution ofy{tQ). Suppose that 

fi^ ^ fi and fi^ ^ fi as e — > 0. (3.22) 

Then, as e ^ 0, 

and ^ u. (3.23) 

(Here stands for the vague convergence.) 



3.2.3 A coupling for walks on different scales. 

The major pitfall of the preceding theorem is that it works only with sequences of 
deterministic speed measures. We want, however, to consider random speed measures 
T^. As we have already remarked, it is not difficult to see that converge to p vaguely in 
distribution. However, it is not enough to make an application of Theorem 13 . 71 p ossible . 
Here the coupling whose existence is stated in Proposition 13.31^ 1) comes into play. It 
allows to replace the convergence in distribution by the almost sure convergence. Then 
it is possible to apply Theorem 13.71 Let us construct this coupling. 
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Consider a two-sided Levy process (or a-stable subordinator) U{x), x G M, f/(0) = 0, 
with stationary and independent increments and cadlag paths defined by 

POO 

^^g-A(c/(x+xo)-(/(xo))] =exp xa / (e"^"" - I)u7'^~°du7 . (3.24) 

Jo ^ 

Let p be the random Lebesgue-Stiehjes measure on M associated to U, p{a, b] = U (b) — 
U{a). It is a known fact that p(dx) = '^jVjS^.^dx), where {xj,Vj) is an inhomogeneous 
Poisson point process with density dx av~^~"dv, which means that p has the same 
distribution as p. 

For each fixed e > 0, we will now define the sequence of i.i.d. random variables such 
that r^'s are functions of U and have the same distribution as tq. Let G : [0, oo) i— > [0, oo) 
be such that 

P(t/(1) > G{x)) = P(ro > x). (3.25) 

It is well defined since U{1) has continuous distribution, it is nondecreasing and right- 
continuous, and hence has nondecreasing right-continuous generalised inverse G^^. 



Lemma 3.8. Let 



G-^ (£-^/°(f/(£(a: + 1)) - U{ex))) . (3.26) 
Then for any e > 0, the are i.i.d. with the same law as tq. 

Proof. By stationarity and independence of increments of U it is sufficient to show 
]P(t^ >t)= P(ro > t). However, 

F{t^ >t)= F{U{e) > e^/"G(t)) (3.27) 

by the definitions of Tq and G. The result then follows from ()3.25|) and the scaling 
invariance of U: U{e) ^^/"[/(l). □ 

Let us now define the random speed measures using the collections {t^} from the 
previous lemma, 

f'{dx) = Y,Cer^^ei{dx). (3.28) 

Finally, using instead of r^, we define the scaling functions S"^ similarly as in ()3.18|) 
and ()3.19|) . The process is then given by = X{t^, S"^), and the construction of 
the coupling from Proposition 13.31 is finished. 



3.2.4 Scaling limit 

Using the three tools introduced above, we can now sketch the proof of Theorem 13.41 
Actually, not many steps remain. 

First, to prove the convergence of to Z it is sufficient to verify the a.s. convergence 
of the speed measures S'^ o to p and then apply Theorem 13.71 The proof of this 
convergence is not difficult, however, slightly lengthy. It can be found in |FINn2j and 
BCn5l. 
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Proposition 3.9. Let and p be defined as above. Then 

S'of'^p and S'ot'^p as e 0, p-a.s. (3.29) 

Finally, to pass from the convergence of X'^ to the convergence of X"^ it is necessary 
to control the scaling functions S^. 

Lemma 3.10. As e —>■ we have 

S%e[e~'y\)-^y, p-a.s., (3.30) 

uniformly on compact intervals. 

Observe that this lemma also implies that the embedded discrete-time random walk 
Y converges, after a renormalisation, to the Brownian motion, independently of the 
value of a. This is valid also if a > and the discrete-time embedded process Y is not 
a simple random walk but a random walk in random environment. 

Since Lemma (3.101 is one of the key parts of the proof of the scaling limit for a ^ 
we prove it here. 

Proof of Lemma \S.l(A We consider only y > 0. The proof for y < is very similar. By 
definition of 5*^ we have eS^{[e~^y\ ) = e X]}lo ^i' where for fixed e the sequence is 
an ergodic sequence of bounded positive random variables. Moreover, ff is independent 
of all fj with j ^ {i — l,i,i-\-l}. The p-a.s. convergence for fixed y is then a consequence 
of the strong law of large numbers for triangular arrays. Note that this law of large 
numbers can be easily proved in our context using the standard methods, because the 
variables f| are bounded and thus their moments of arbitrary large degree are finite. 
The uniform convergence on compact intervals is easy to prove using the fact that is 
increasing and the identity function is continuous. □ 



3.3 Aging results 

The aging results for the one- dimensional BTM follow essentially from Theorem 13.41 To 
control the two-time functions R and it is only necessary to extend its validity to the 
joint distribution of (X^(1),X^(1 + 6)) at two fixed times, which is not difficult. This 
extension then yields the following aging result. 

Theorem 3.11 (Aging in the one-dimensional BTM). For any a E (0, 1), > 

and a G [0, 1] there exist aging functions Ri{6), R'^{9) such that 

lim R{t^,t^ + et^) = lim m>[x{{i + d)t^) = x{t^)\T] = R^{e), 



lim R'^it^, t^ + et^) = lim E y [P(X((1 + 9)t^) = i\t, X{t^))]^ = R^O). ^^'^^^ 



Moreover, Ri{6) and R'^{6) can be expressed using the analogous quantities defined using 
the singular diffusion Z : 

R,{9) =EF[Z{l + 9) = Z{l)\p], 

Ri{e) = E^[P(Z(l + e) = x\p, Z{1))]\ (3-32) 
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Remark. 1. This result is contained in |FINn2j for a = and in (BCOSj for a > 0. 

2. Let us emphasise that the functions Ri{6), R'^{9) do not depend on the parameter 
a, since the diffusion Z{t) and the measure p do not depend on it. This is the result of 
the compensation of shorter visits of deep traps by the attraction to them. 

3. It should be also underlined that only averaged functions are considered in the 
theorem. For a fixed realisation of r there is no limit of + Otu,'-, t). From the 
proof of Theorem 13. lit it is however not difficult to derive the following weaker result 
(see Theorem 1.3 in jCerOGj ). 

Theorem 3.12 (Quenched aging on Z, in distribution). As — > oo, the distri- 
bution of R(tw,tw + Otnj] t) converges weakly to the distribution ofF[Z{l + 9) = Z{l)\p]. 

3.4 Subaging results 

In the case of the two-time function 11 much shorter times t should be considered, t <^ t^, 
i.e. subaging takes place: 

Theorem 3.13. For any a G (0, 1), 6* > and a G [0, 1] there exist an aging function 
ni,a(6') such that 

\imu{t,t + fa{t,e)) = \imEF[x{t') = x{t)yt' e [t,t + fa{t,e)]\T] =ni,,(^^), (3.33) 

where the function fa is given by 

fa{t,9) = 9t^^'-^'^L{tY-\ (3.34) 

Here we use 7 to denote the subaging exponent, 7 = (1 + a)~^ , and L(t) is a slowly 
varying function that is determined only by the distribution of tq. The function lii^a{0) 
can again be written using the singular diffusion Z , 

POO 

ni,a(^) = / glieu'^~')dF{u), (3.35) 
Jo 

where F{u) = EP[p({Z(l)}) < u\p], and where ga{\) is the Laplace transform of the 
random variable VaTQ, 

(7,(A)=E(e-"'^"^o). (3.36) 

If a = 0, ()3.35p can be written as 

POO 

ni,o(^) = / e-^/MF(M). (3.37) 

Remark. 1. As can be seen, in this case the function Ili a{d) depends on a. This is not 
surprising since the compensation by attraction has no influence here and the jump rates 
clearly depend on a. 

2. Of course, an analogous result to Theorem 13.121 holds for quenched subaging in 
distribution. 
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In the RHT case, a = 0, the proof of Theorem 13 .131 is straightforward. It follows from 
the convergence of the distribution of 'rx^(i) to the distribution of p({Z(l)}) as stated in 
Theorem 13.41 This then implies the convergence of the distribution of Tx{t^) I {t1jL{tw)) ■ 

The proof in the case a > is more complicated. It essentially involves a control 
of the distribution of the depth of those traps that are nearest neighbours of the traps 
where X(t^) is with a large probability. It turns out that this distribution converges to 
the the distribution of tq, i.e. there is nothing special on the neighbours of deep traps. 

The behaviour of the two-point functions n(t„,, t + tw) and R{tw,t + tw) is not difficult 
to understand and guess. We give here a heuristic explanation for these results, first in 
the case a = 0. After the first n jumps the process typically visits 0(n^/^) sites. The 
deepest trap that it finds during n jumps has therefore a depth of order 0(n^/^°), which 
is the order of the maximum of n^/^ heavy-tailed random variables as can be verified 
from Assumption 12.31 This trap is typically visited 0(n^/^) times. Since the depths are 
in the domain of attraction of an a-stable law with a < 1, the time needed for n jumps 
is essentially determined by the time spent in the deepest trap. This time is therefore 
Q(^^{i+a)/2a^_ Inverting this expression we get that the process visits typically 0{t°''^) 
sites before time t. The deepest traps it finds during this time have a depth of order 
f^. Moreover, the process is located in one of these deep traps at time t. From this, one 
sees that the main contribution to quantity R{tw^tw + 1) comes from the trajectories of 
X that, between times t^, and tw + 1, leave the original site X(t^) a number of times of 
order , and then return to it. Each visit of the original site lasts an amount of time 
of order t^,) which is the time scale on which 11 ages. 

For a > 0, such heuristics is not directly accessible. However, the Theorem 13 . 41 vields 
that Tx{tn,) is of the same order, 0(^2,), as in the RHT case. Each visit of X(t^) lasts a 
shorter time, 0(tX''^~"^), as follows from the definition of the process and the fact that 
the depths of the neighbours of X{tyj) are 0(1). On the other hand, the process makes 
more excursions from X{tw), their number being 0(tw"^"''). 

3.5 Behaviour of the aging functions on different time scales 

Having found two interesting time scales t = O(t^) and t = 0{t1j) in the model^ one may 
ask if there are other interesting time scales for the functions H and R. This question 
was raised by Bouchaud and Bertin in |BBn3j . The negative answer was given in |Cer06j : 

Theorem 3.14 (Behaviour of H on different scales), (a) Short time scales. Let 

f(t) he an increasing function satisfying > f(t) > for allt large and for some 
J > K, > fj, > 0. Then 




(3.38) 



with < i^i < oo. 



^We suppose here that a = and Assumption 12 . 31 holds . that is L{t) 1. 
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(b) Long time scales. Let g{t) be such that = o{g{t)). Then 



^lim (^)"n(t,t + ^7(t)) =^^2, (3.39) 

with Q < K2 < oo. 
(c) Behaviour of ni^o(^)- The function 11(6') satisfies 

limr-^(l -nio(^)) = i^i, (3.40) 

lim rni,o(^) = K2. (3.41) 

Remark. We emphasise that the constant Ki occur both in ()3.38p and ()3.40|) . That 
means that the behaviour of Hi 0(6*) at 6* ~ gives also the behaviour of Il(t,t + fit)) 
for f{t) -C f^. An analogous remark applies for time scales f{t) ^ f^. Both constants 
Ki and K2 can again be expressed using the F.I.N, diffusion. The formulae , and a lso 
similar (slightly weaker) results for the two-time function R can be found in (CerOB . 



We give again a heuristic description of this result. As we know, at time t the process 
X is typically in a trap of depth Olf), it needs a time of the same order to jump out. In 
Theorem 13. 14r a) we look at 1 — Il{t,t + f{t)) with f{t) -C f^, that is at the probability 
that a jump occurs in a time much shorter than f . There are essentially two possible 
extreme strategies which lead to such an event: 

1. Tx(t) has the typical order but the jump occurs in an exceptionally short time. 
X{t) is in an atypically shallow trap and stays there a typical time. 



In |Cer06j it is proved that the second strategy dominates. Therefore, one has to study 



the probability of being in a very shallow trap or, equivalently, to describe the tail of 
^[Tx{t)/t'^ < u] for M ~ 0. To control this tail the proof use the fact that although the 
BTM never reaches equilibrium in a finite time in infinite volume, it is nearby equilibrium 
if one observes only traps that are much shallower than the typical depth on intervals 
that are small with respect to the typical size of X{t). This puts on a rigorous basis, at 
least in dimension one, the concept of local equilibrium that was introduced in the physics 
literature by [RMBnOj . The concept does not give the right predictions for the values 
of the hmiting functions -Ri(^) and Hi 0(6') but it is useful to describe their asymptotic 
behaviour. A very similar heuristics applies also for time scales f(t) S> f^. 



3.6 References 

In the RHT case, a = 0, this model has first time been studied by Pontes, Isopi 
and Newman in |FIN99j . It was used there as a tool to control the behaviour of the 
voter model with random rates. It was proved there that the process is sub- diffusive, 
i.e. E[X(t)/-\/t] as t — 00, and that the dynamics localises in the sense that (for 
a.e. r) 

supP[X(t) = x|t] 74 ast^oo. (3.42) 
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That means that there will be always a site (dependent on time t) where the process 
X can be found with a non-negligible probability. This localisation occurs only in one- 
dimensional BTM and is at the heart of the majority of the differences between Z and 

other graphs. 

Th e scaling limit has been established for the case a = in |KINn2| and bv ICerOal 
BCn5 for a 7^ 0. Following |BBn3j . the results of Section ITHl have been given by |Cer06j . 



For aging in another interesting one-dimensional dynamics, i.e. Sinai's Random 
Walk, see |DGZOH IBFOGj . Let us mention two interesting open questions on this one- 
dimensional trap model: 

(a) What is the behaviour of the edge of the spectrum for the generator of the 
dynamics. This might be close to, but easier than the same question solved for Sinai's 
Random Walk by |HKn6j . 

(b) What is the influence of a drift in the BTM? Monthus |Monn4j gives a very 
interesting picture based on renormalisation arguments. 



4 The trap model in dimension larger than one 

After resolving the BTM on Z, the next natural step is to study the Bouchaud model 
on the (i- dimensional lattice, BTM(Z'^, r, 0), d > 1. Observe that we set a = 0, that 
means that only the RHT dynamics is considered. In this section we always assume that 
Assumption 12.31 holds. 



4.1 The fractional-kinetics process 

As for the one- dimensional model, we first identify a scaling li mit of th e BTM on Z"^. 
The result of this section is contained in the forthcoming paper jBCOBbj . We will, from 
now on, use frequently the theory of Levy processes and subordinators. A very short 
summary of this theory can be found in Appendix El 

Let us first define the process that appears as the scaling limit. 

Definition 4.1 (Fractional kinetics). Let B^it) be the standard (i-dimensional Brow- 
nian motion started at and let V be the a-stable subordinator given by its Laplace 
transform E[e-^^W] = e"*^". Let T(s) = inf{t : V{t) > s} be the inverse of V{t). We 
define the fractional-kinetics process by 

= 5,(r(s)). (4.1) 

We list here without proofs several properties of the process 

Proposition 4.2. 1. The stable subordinator V is strictly increasing, therefore its in- 
verse T and also the process are continuous. 

2. The name of the process is due to the following fact. Let p{t, ■) be the probability 
density of'^dit). Then p is a solution of the fractional kinetic equation, 

^p(t,x) = -Ap(t,x) + 6(0)—'——. (4.2) 
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Here, the fractional derivative d°'p{t, x) / dt" is the inverse Laplace transform of s"'p{s, x), 
where p{s, x) = e~^^p(t, x) dt is the usual Laplace transform. The equation ()4.2|) 
should be understood in the weak sense, i.e. it holds after the integration against smooth 
test functions. 

3. The process is not Markov, as can be seen easily from the previous point. 

4. The fixed-time distribution of'^d is the Mittag-Leffler distribution, 

E(e*«-*''W) =E,(-|e|V), (4.3) 

where Ea{z) = Z]m=o -^"^/^(l + ^0). 

5. The process \E'd is self-similar: 

m-^d) '= A-°/'$(At). (4.4) 

The process '^d is well known in the physics literature (see |Zasn2j for a broad survey 
and earlier references). It is the scaling limit of a very classical object, a Continuous 
Time Random Walk (CTRW) introduced by |MW65j . More precisely consider a simple 
random walk Y on and a sequence of positive i.i.d. random variables {sj : i G N} 
with the distribution in the domain of attraction of an a-stable law. Define the CTRW 
U{t) by 

fc-l k 

U{t) =Y{k) iite[J2^^^Yl ^0 • ^^-^^ 

i=l i=l 

It is proved in .MSn4j that there is a constant C (depending only on the distribution of 
Si) such that 

Cn-'^'^Uitn) ^rf(t). (4.6) 



4.2 Scaling limit 

Observe that the Bouchaud model (for a = 0) can be expressed as a time change of the 
simple random walk. The time-change process is crucial for us: 

Definition 4.3. Let 5(0) = and let S{k), A; e N, be the time of the fc*'' jump of X. For 
s G M we define S{s) = S{[s\). We call S{s) the clock process. Obviously, X(t) = Y{k) 
for all S{k) <t<slk + l). 

The following result shows that the limit of the d-dimensional Bouchaud model and 
its clock process on Z'^ {d > 2) is trivial, in the sense that it is identical with the scaling 
limit of the much simpler ( "completely annealed" ) dynamics of the CTRW. 

Theorem 4.4 (Scaling limit of BTM on Z^). Let 

|n"/2, tfd>3. 
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Then for all d > 2 and for a.e. r 



fin) 



n- 



■oo 



^ ^d(t) and 



S{C,{a)-'f{n)H) 



71- 



■OO 



V{s). 



(4.8) 



n 



weakly in the Skorokhod topology on D{[0,T],'R'^) (the space of cadlag functions from 
[0, T] to M.'^). If Gd{0) denotes Green's function of the d- dimensional random walk at 0, 
then 



The main ideas of the proof of this theorem will be explained in Section 14.41 At 
this place, let us only compare the fractional-kinetics process with the F.I.N, dif- 
fusion Z. Both these processes are defined as a time change of the Brownian motion 
Bd{t). The clock processes however differ considerably. For c/ = 1, the clock equals 
0(t) = / £(t, x)p(dx), where p is the random speed measure obtained as the scaling limit 
of the environment. Moreover, since i is the local time of the Brownian motion Bi, the 
processes Bi and are dependent. For c? > 2, the Brownian motion Bd and the clock 
process, i.e. the stable subordinator V , are independent. The asymptotic independence 
of the clock process S and the location Y of the BTM is a very remarkable feature dis- 
tinguishing d > 2 and d = 1. It explains the "triviality" of the scaling limit in dimension 
d > 2, but is, by no means, trivial matter to prove. We will come back to an intuitive ex- 
planation of the independence in Section 15. HI Note also that nothing like a scaling limit 
of the random environment appears in the definition of "^d, moreover, the convergence 
holds T-a.s. The absence of the scaling limit of the environment in the definition of "^d 
transforms into the non-Markovianity of "^d- Note however that it is considerably easier 
to control the behaviour of "^d than of Z even if "^d is not Markov: many quantities 
related to \E'rf can be computed explicitly, as can be seen from Proposition 14.21 

4.3 Aging results 

The following two theorems describe the aging behaviour of the two-time functions R 
and n. 

Theorem 4.5 (Quenched aging on Z'^). For all a G (0, 1) and d > 2 there exists a 
deterministic function R{0) independent of d (but dependent on a) such that for P-a.e. 
realisation of the random environment r 



The function R{-) can be written explicitly: Let As\a{u) be the distribution function of 
the generalised arcsine law with parameter a, 




(4.9) 



lim R{tu,, tyj -\- Otyj] T 



)=R{9). 



(4.10) 




(4.11) 



Then R{9) = AsU(l/l + ^). 
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Theorem 4.6 (Quenched (sub-)aging on l/). For all a E (0,1) and d > 2 there 
exists a deterministic function lld{0) such that for P-a.e. realisation of the random en- 
vironment T 

lim u{t^,t^ + ef{t^y,T) = iiaie), (4.12) 

with 

ifd = 2, 



/(tj = <( "'^ (4.13) 
tw, if d > 3. 

The function Ud does depend on d. Moreover, for all 6 > it satisfies 

lim Ud{e) = R{e). (4.14) 

d— >oo 



Remark. 1. Both theorems are proved in BCMOGj for d = 2 and in Cern3j for d > 3. 



The proofs are relatively technical and exceed the scope of these notes. The main ideas 
however do not use specific properties of the integer lattice Z'^ and can be generalised to 
different graphs. These ideas will be explained in Section The other important part 
of the proof, that is the coarse-graining of the trajectory of the process, is explained in 
the next subsection. 

2. The function 11^ can also be explicitly calculated but the formula is tedious 



BCM06 Cern3 



3. We will see later that the function R is closely related to the arcsine law for Levy 
processes. As we will also see, the same function appears as the limit in the case of the 
BTM on a large complete graph. Therefore, the mean-field dimension of the BTM for 
the two-time function R is d = 2. 

4. Both presented results are quenched, i.e. they hold for a.e. r. This should be 
compared with the results for the d = 1 case f Theorems 13. lit 13.121 and 13.13^ where only 
averaged aging holds. The reason for this difference is explained in Section 1331 below. It 
is, of course, trivial to get averaged results from Theorems 14.51 and 14.61 the dominated 
convergence theorem yields 

lim R{t^,t^ + et^) = R{e), 

, , (4.15) 

lim u{t^,t^ + ef{t^)) =u{e). 

5. The fact that in > 3 the time scale f{tw) = t^, is the same for both functions 
R and 11 is a consequence of the transience of the simple random walk. The traps are 
visited only a finite number of times, so that different time scales cannot appear for R 
and n. 

4.4 The coarse-graining procedure 

We would like to describe here the coarse-graining procedure which was introduced in 



BCM06 , and which is the main tool in proving Theorems I4.4H4.6I Even this short 



sketch of the procedure might be considered technical and can be skipped on a first 
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reading. The reader can also decide to return here after being acquainted with the 
general ideas of Sections 15.21 and I5.31 We will deal here only with the convergence of the 
clock process to an a-stable subordinator. For the sake of concreteness we set d = 2, 
however the same arguments apply also for d > 2. The discussion in this section is valid 
for a.e. realisation of the random environment r. 

The random environment is heavy tailed. Therefore, as for d = 1, the behaviour of 
the clock process S{k) is determined by the time spent in the deepest traps that the 
process visits during the first k steps. It is thus necessary to find the depth scale of these 
traps and then study how these traps contribute to the clock process. 



The coarse-graining procedure of jBCMn6j studies the process Y (resp. X) only before 



the exit from a large disk I]){n), n G N, of area m2^n^~°' centred at the origin. We denote 
by R{n) the radius of this disk, i?„ = ^/7r~^2^n^~". The random walk Y makes 0{R^) 
steps in D(n) and visits 0(i?^/logi?^) different traps, as is well known. Therefore, the 
deepest traps that it visits have a depth of order 0((i?^/ log i?^)^/") = 0(2"'/"n~^). We 
therefore approximate the clock process by the time that X spends in the set 

T^' = {xe D(n) : £2"/"^-^ < < M2"/"n-^}, (4.16) 

where e is a small and M a large constant. We call the traps in this set the deep traps. 
It must be proved that this approxim ation is correct. We do not want, however, to deal 
with this problem here (see (BCMOGj ). 

After giving the proper depth scale, we can study how the deep traps are visited. 
This is where the coarse-graining is crucial. We cut the trajectory of the process Y into 
short pieces. Every such piece of the trajectory ends when Y exits for the first time the 
disk of area 2'^rC around its initial point. At this moment a new piece starts. Clearly, 
we should take 7 < 1 — a. Formally, we set j'q = 0, and then we define recursively 



= min {A; > jU ■ dist(F(A;), Y{jU)) > vV^2^}. (4.17) 

We use to denote the starting points of the pieces of the trajectory, = Y{jf). It 
can be seen easily that the number of pieces of the trajectory before the exit from D(n) 
is of order 0{n^~°'~'^). 

We look at the time that the walk spends in the deep traps during one piece of 
trajectory: we define the score of the piece i by 

Ji + l ^ 

E efcry(,)ll{r(fc)GTf}. (4.18) 

To study the scores it is convenient to introduce another family of random variables 
indexed by x G D(n). We set the distribution of to be the same as the distribution of 
conditioned on the fact that the i*'^ piece of the trajectory starts at x, i.e. conditioned 
on = x. The main technical piece of the proof is to show that the law of does not 
depend on x (with a small exceptional set): 
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Lemma 4.7. Fix k > large enough and define 



£(n) = {x e B{n) : dist(a;,T/^) > Vn-^^hF^}. (4.19) 
Then for F-a.e. random environment t, uniformly for x G £{n) 



Here 



a 1 



F(A) = F(A; M, a) = /C ( pf - / " ) (4.21) 



with K,' = 7r-Mog2 and K. = (log2)-i. 

This lemma is a consequence of the following four facts, whose proofs are based on 
classical sharp estimates on the Green's function for the simple random walk on Z*^ and 
on certain "homogeneity" properties of the random environment. These proofs can be 



found in BCMOB 



1. is equal to with probability l-/Cpfn"+T-^(l + o(l)), where = e-^-M"". 
That means that typically no deep trap is visited in a piece of the trajectory. Further, 
it implies that only a finite number of pieces has a non-zero score before the exit from 
D(n). 

2. With probability )Cp^n°'^'^^^{l + o(l)) the random walk visits (many times) only 
one deep trap during one piece, call it y. The probability that two or more deep traps 
are visited during one piece is 0{n^^°'~^'^~^^), therefore, with overwhelming probability, 
this event does not occur before the exit from D(n). 

3. In the case when one deep trap y is visited, the distribution of its normalised 
depth 2~"'/"nry converges to the distribution on [e, M] with the density p{u) proportional 
to u~°'-^. 

4. The number of visits to y is geometrically distributed with mean }C'n{l + o(l)). 
Therefore, conditionally on hitting y the score has an exponential distribution with mean 
KlnTy which is of order 0(2"/°). 

Using Lemma ITTj it is not difficult to check that the scores are asymptotically i.i.d. 
Actually, to transfer the uniform convergence of distributions of s" into the asymptotic 
i.i.d. property of it is sufficient to check that with an overwhelming probability all 
pieces of trajectory before the exit of ©(n) do start in £{n). 

The asymptotic i.i.d. property then yields the convergence of the normalised sum of 
scores, 2~"/" ^^"o ^ "^^5 to a Levy process. The Levy measure of this process can be 
computed from Lemma 14.71 Using the knowledge of the Levy measure, one can then 
prove that for any T > it is possible to choose m large enough such that as e — 
and M — > oo the distribution of this Levy process on [0, T] approaches the distribution 
of an a-stable subordinator on [0,T]. The convergence of the clock process to the same 
subordinator then follows since the sum of scores is a good approximation of the clock 
process. 
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4.5 References 



]. In the case a = 0, 
. This proof is based 



The trap model in dimension > 2 is studied in IMB961 IRMBO 
a mathematical proof of aging has been given in jCerOSl lBCMn6| 
on the c oarse-graining of the trajectories of the BTM sketched in the last section. In 
fBC06a we establish the fractional-kinetics scaling limit based on the arguments of 
^BCM06j . The case where a 7^ is discussed in |RMB01j . but is still an open problem. 

Let us mention that the clock process introduced in Definition 14.31 is close to the 
problem of Random Walk in Random Scenery (RWRS), except for the extra randomisa- 
tion due to exponential waiting times. However, here the tails of the scenery distribution 
are heavier than in the recent works on RWRS. 

In dimension d = 1, RWRS with heavy tails have been studied by Kesten and Spitzer 
jESZi- The F.I.N, diffusion could indeed be seen as a Brownian motion time-changed 
by a (dependent) Kesten-Spitzer clock. 



5 The arcsine law as a universal aging scheme 

In this section we explain a general strategy that can be used to prove aging of the 
functions 11 and R in the BTM on many different graphs (including Z"^ for d > 2, tori 
in Z*^, large complete graphs, and high-dimensional hypercubes). When this strategy 
can be used, the behaviour of R{tw, + 9tw', t) for large can be expressed using the 
distribution function of the generalised arcsine distribution with parameter a, AsIq,(-). 

More formally, we will consider in this section a sequence of the Bouchaud trap 
models {BTM(G„, t„, 0) : n G N}. We will prove that for properly chosen time scales 
twin) 

lim RnitM, (1 + 0)tM; t) = AsU(l/l + 6), (5.1) 

n— >oo 

where 

Rnit^,t^ + t;T) ■.= ¥[XM=Xnitw + t)\T]. (5.2) 

The possible time scales tu,(n) will depend on the graphs Gn and the laws of the depths 
r„. We will always try to prove aging results on the widest possible range of time scales. 

We will consider only the case of RHT dynamics, i.e. a = 0. So that, the embedded 
discrete-time process Yn is the simple random walk on G„, and the continuous time 
Markov chain X„ is a time change of Y. Our strategy relies on an approximation 
of the clock process (see Definition 14. 3j) . more precisely of its rescaling, by an a-stable 
subordinator. We will then show that the event Xnitw) = Xn{{l + 0)tw) is approximated 
by the event that the subordinator jumps over the interval [1, 1 + ^]. The classical arcsine 
law for Levy processes (see Proposition I A. 41 in the Appendix) will then imply the aging 
resuh ira . 

5.1 Aging on large complete graphs 

As a warm-up, and in order to explain our strategy, we give here a complete proof of 
aging on a large complete graph in this sub-section. The advantage of this graph is that 
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the embedded simple random walk is particularly simple. On the other hand, almost 
all effects of the trapping landscape are already present. The method of the proof that 
we use here is probably not the simplest one. Its main ideas can however be adapted 
to more complex graphs; the structure of the proof stays the same, but a relatively fine 
control of the simple random walk on these graphs is then required. 

Note that the complete graph is the graph that was proposed in the original paper 
of Bouchaud |Bou92j . The aging on this graph is proved using renewal arguments in 
BD95j . see also the introduction to |BBG03bj . Another proof of aging and much more 



can be found in |BF05j where eigenvalues and eigenvectors of the generator of X are 
very carefully analysed. 

Let Gn be a complete graph with n vertices, Gn = (V„, £n), where V„ = {1, . . . ,n} 
and £n = {x,y),x,y G V„. Note that we include loops {x,x) into the graph, so that 
jumps from x to a; are possible. This makes the embedded random walk of the RHT 
dynamics extremely simple: the positions Yn{i), i E N, are i.i.d. uniform random vari- 
ables on V„. We also suppose that the starting position of the process, X„(0) = Yn{0), 
is uniformly distributed on Vn- Finally, we assume that the trapping landscape is given 
by an i.i.d. sequence {t^ : x eN} independent of n such that Assumption 12. HI holds. 

As the process Xn is a time change of the i.i.d. sequence Yn{i), we need to study the 
clock process Sn(t). Observe that Sn{k) for G N is the time of k^^ jump of Xn, 

k-l 

Sn{k) = '^TY„{i)ei, (5.3) 

where {cj} is an i.i.d. sequence of exponential random variables with mean one. 

Since we have included the loops into the graph, we redefine the two-time function 
n slightly: 

n„(t^, + t; r) := F[{Sn{k) :keN}n [t^, + t] = 0|t] . (5.4) 

Note that this definition differs from the original one only if a jump from X{tw) to 
X{tw) occurs. The probability of this event is and is thus negligible for large n. 
Similarly, for the complete graph the function i?„ differs very little from the function 
n„, a correction is again of order 1/n. 

We prove aging on time scales smaller than n^/". The scale n^/" appears because the 
deepest trap in Vn has a depth of this order as can be easily verified from Assumption 12 . 31 
Therefore, at time scales shorter than n^^" the process has not enough time to reach the 
equilibrium, and aging can be observed: 

Theorem 5.1 (Quenched aging on the complete graph). Let < k, < 1/a and let 

tw = tw{n) = . Then for a.e. random environment r 

lim Iin{tu,{n), (1 + e)tM\ r) = AsU(l/l + 9). (5.5) 

n— >oo 

Moreover, the clock process converges to the stable subordinator V with the Levy measure 
ar(l + a)u-"-MM, 



n 



K,/a 



27 



Remark. 1. A similar result holds on the shortest possible time scale k = 0: if = 
tn^ = t, then 

lim lim Unit^, (1 + 6%; r) = AsU(l/l + 9). (5.7) 

This is actually proved in jBD95j and |BFn5j . 

2. For the longest possible time scales tw{n) = tn^^" a double limiting procedure is 
also necessary. Essentially the same arguments that are used to prove Theorem 15 . 1 1 vield 
that n„(t^(ri), (1 + 6)t^{n); r) converges to AsIq,(1/1 + 6) in probability, that is for any 
e > 



lim lim P 

t— >0 n— »oo 



n(tn^/",(l + ^)tn^/";T) - AsU(l/l + 0) > £ = 0. (5. 



3. No scaling limit for X exists since there is no such thing as a scaling limit for the 
complete graph. 

To prove Theorem 15.11 we consider a rescaling of the time change process, 

Sn{t) = ^Snik) for t e [un{k), u^k + 1)), (5.9) 

where M„(fc) = Yli=i {e^, i G N} is another sequence of mean one i.i.d. expo- 

nential random variables. We introduce this sequence in order to randomise the times 
of jumps of Sn- Thanks to this randomisation Sn is the Levy process. The standard 
rescaling of Sn, n~'^Sn{tn'^") does not have this property. Sn is compound Poisson pro- 
cess. The intensity of its jumps is n°"^. Every jump has the same distribution as n~'^Txei, 
where x is uniformly distributed in V„. Therefore, the Levy measure fin of Sn is 

fin{du) = n°"- V — du = ^""+''-1 V du. (5.10) 

xeVn xeVn 

It follows from definitions of n„ and iS„ that 

^„(n^ (1 + e)n^; t) = P[{5„(t) : t E R} D [1,1 + 6] = (5.11) 

that is the probability that Sn jumps over the interval [1, 1 + 6*]. 

The idea of the proof is the following. We write Sn as a sum of three independent 
Levy processes, Sn,M = S"^^, S^^ and Sf^ = S^q with the Levy measures /i„,M = I^^mi 
l^n,e and /i^ = /i^^o) where 

<,(d«) = n--+^-' e T,^(n)}le-"'^"/- du, (5.12) 

and 

T,V) = {a:GV„:^G[a,6)} (5.13) 

That is we divide traps into three categories: (a) the very deep traps, x G TM{n), 
i.e. Tx > Mn'^, (b) the deep traps, x G T^'^{n), i.e. Tx/n'^ G [e,M), (c) the shallow traps. 
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X G T"^, i.e. Tx < en'^. We consider the contributions of these different categories to the 
clock process separately. 

We show that the contribution of the deep traps, i.e. the process iSf^, is well ap- 
proximated by an a-stable subordinator, at least if e is small and M large enough. The 
proof of this fact relies on the weak convergence of the Levy measures /if,. We prove 
this convergence only on the interval [0,T] where T is chosen such that iSf^ is larger 
than (1 + 6) with a large probability. 

Further we prove that the very deep and shallow traps can be almost neglected. 
More exactly, we prove that Sn^niT) = with a large probability, that is Sn,M does 
not jump before T, and that Sf^{T) can be made small by choosing e small enough. 
These fact together will imply that Sn is well approximated on [0, T] by an a-stable 
subordinator, and the claim of the theorem follows from ()5.1H1 and the arcsine law for 
stable subordinators, more precisely from CoroUarv IA.5I 



5.1.1 Deep traps 

The following proposition describes the contribution of the deep traps to the time change, 
that is the process 5^^^. 

Proposition 5.2. The Levy measures /xf^ of S^^ converge weakly as n ^ oo to the 

measure /xf given by 

fir{du)=l ^e-^/^dzdu. (5.14) 



This proposition has an important corollary that states that S^\. can be well approx- 
imated for large n by an a-stable subordinator. To this end, let Z^^ be a subordinator 
with the Levy measure 

/if (d«) = e-"/^ dz du+ e-"/^ dz du (5.15) 

Jo ^ Jm z 

independent of all already introduced random variables. 

Corollary 5.3. The process iSf^ + ^f converge weakly in the Skorokhod topology on 
i5([0,T],]R) to an a-stable subordinator. Moreover, for any T and 6 > it is possible to 
choose e small and M large enough such that 

P[Zf (T) >6]<6. (5.16) 

Proof of Corollaru \5.^ The first claim is consequence of the convergence of Levy mea- 
sures (Proposition 15. 2j) together with Lemma IA.21 The fact that the limit is a stable 
subordinator follows from 

/if (du) + /if (dw) = ar(l + a)^-"-^ du. (5.17) 

Finally, to prove ()5.16|) observe that for all A > 

^e-^^riT)^^ = g-Tj- (i-e--)Af (dx) e-^OM^ ^^^^g^ 

that is the law of Z^\T) converges weakly to the Dirac mass at 0. □ 
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To prove Proposition 15.21 we first show that for a.e. r there is the right number of 
deep traps with depths approximately un'^ in V„, m G [e,M). This follows from a "law- 
of-large-number-type" argument using the fact that we have a large number of traps 
with such depth. This property fails to be true if k > 1/a (which explains why we 
cannot get a.s. result for ty^ = cn^^" and why we need the double limit procedure). 

Lemma 5.4. Let < k < 1/a. Then there exists a function h{n), h{n) ^ as n ^ oo, 
such that 



lim sup 

ue[e,M)n{Z/hin)} 



u+h{n) 



T-a.s. 



(5.19) 



Proof. Let g{u) be defined by F[tx > u] = u~°'{l+g{u)). It follows from Assumption 12.31 
that lim„^oofi'(w) = 0. Take h{n) such that h{n) — » 0, h{n) > (logn)~^, h{n) ^ 
sup„g[g^^,j) (7 (■un''). Then, as n — > oo, 



P[r^ e (u + h{n))n'')] = an"^-"/i(n)n-"'^(l + o(l)). 
Indeed, by definition of g, the left-hand side is equal to 



(5.20) 



u ) 



g[un )-g 



{{u^h{n))n^)(\^^^ (5.21) 



The last expression is equal to the right-hand side of ()5.20j) as follows from the definitions 
of g and h. 

Fix M G [e, M] and £ > 0. Using twice the exponential Chebyshev inequality together 
with (IH:^ we get 



P 



|Tn"+^("^(n)| 



I'l — e, 1 + e)au 



-o-l 



< Cexp{-c(e)M-"-ini-""/i(n)}. (5.22) 



Summing over u G [e, M) fl {Z//i(n)}, the claim of the lemma then follows using the 
Borel-Cantelli Lemma and the fact that h(n) > (logn)~^. □ 

Proof of Proposition [5~^ By definition of /x^^^. 



n 



n 



xdTM (n) 
{M-e)/h(n) 



i=0 rp{i + l)h{n) 

ih(n) 



{M-s)/h{n) 



- E 



IT. 



(i+l)h(n), 



ih{n) 



n] 



i=0 



n'^ih{n] 

,, 1 
ih{n) 



dM(l + o(l)) 



(5.23) 



-u/ih(n) 



du(l + o(l)). 
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By the previous lemma, for large n, r-a.s 

{M-e)/h(n) -u/ih{n) 

/x5(dn) = <^mn))-^-''---—h{n) d«(l + o(l)), (5.24) 

which is the Riemann sum of the integral in ()5.14jl with the mesh size h{n). Since 
h{n) as n ^ oo, the proof is finished. □ 



5.1.2 Shallow traps 

We prove that the contribution of the shallow traps at any final instant T, S^{T), is 
small if e is chosen small enough. 

Lemma 5.5. There exists a large constant K such that for allT > and n large 

E[5^(T)|r] < fsTTe^"", r-a.s. (5.25) 

Proof. Let ^ be such that T G [un{0,Un{i + 1)) (see ()5.9j) ). It is easy to see that ^ has 
Poisson distribution, E[,^] = Tn'^". By definition of S'^ 



^niT) = ^£e.ry„(,)]l{F„(0 G T^{n)}. (5.26) 

1=0 

To bound its expected value we divide T'^{n) into slices T^2-^^^ ^ i eN. We define 

A^{z) := |r : E [-^ ^ e,ry„(,)]l{F„(0 G T,f_7'''(n)} |t] < irV-°2^("-^)}. (5.27) 



1=0 

We show that there is K' such that a.s. 



An := |^y4„(i) occurs for all but finitely many n. (5.2^ 



i=0 



Lemma (5.51 is then a direct consequence of this claim. 

To show ()5.28|1 we first estimate the probability that a fixed site, say 1, is in T^^_l"^ , 

Pn,i := P[l G T|2_7+'] < C£-°2^"n-^" (5.29) 
as follows from Assumption 12.31 For any x G V„, 



El ' 

' i=0 



[5^e,]l{>;.(0 = x}] = Tn^^-\ (5.30) 
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therefore 



P 



i=0 



< P[ ^ W"-i£2-^+i]l{x e T^2-'^\n)} > i^'ei-°2*("-i)l . (5.31) 



xev„ 



By Chebyshev inequality this is bounded from above by 

exp{-Ani-'^"cir'£-°2^"} JJ (p„,i(e^ - 1) + 1). 



(5.32) 



xeVn 



Using 1 + X < and ()5.29p , it is easy to see that for K' large enough this is smaller 
than exp{— c?T,^~'^"2*"£:~"}. Therefore, summing over i eN, 



(5.33) 



An application of the Borel-Cantelli lemma finishes the proof of ()5.28|) and therefore of 
Lemma 15.51 □ 

5.1.3 Very deep traps 

We show that, with a large probability, none of the very deep traps is visited during the 
first Tn'^" steps: 



Lemma 5.6. If < k, < l/a, then for all n large enough, T-a.s. 

P[5„,Af(T) ^ 0|r] < cTM-" 



(5.34) 



Proof. Define H = mf{k G No : Yn{k) G TM{n)}. By definition of Sn,M, the claim of the 
lemma is equivalent to F[H < C,\t] < cTM~"', with the same ^ as in the previous proof. 
Since 

F[i/ < ^|t] < F[H < 2W"|r] + P[^ > 2W], (5.35) 

and the second term on the right-hand side decreases exponentially in n, it is sufficient 
to bound the first term. Using the fact that r-a.s. \TM{n)\ < CM~°'n^~'^" (which can 
be proved similarly as Lemma f5.4|l . we have 



F[H < 2W"|t] = P IJ {Yn{k) G T^^{n)} 



k=0 



< 2Tn'' 



T^Hn) 



n 



This finishes the proof. 



< cM-"T. (5.36) 
□ 
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5.1.4 Proof of Theorem 15.11 



We can now finish the proof of aging on a large complete graph. All claims here are 
valid T-a.s. 

First, for 5 > we fix T > such that (for all n large enough) 



ns'JAT) + z^\T)<2 + e\T]<-^, 



(5.37) 



which is possible due to Corollary 15. 31 Further, we use this corollary. Lemmas 15 . 51 and l5 . 61 
to fix e and M such that 



P[5„,m(T) ^ 0|t] < ^, 



P 



s->i 



< - and P 
~ 4 



5 



6 

< -. 



2J - 4 



(5.38) 



Let B{n) be the intersection of all events from the two previous displays. It follows 
that P[i?(n)^] < 6. Let E{n) be the event whose probability we are trying to estimate 
(see dSini)): 

E{n) := {{Sn{t) ■.teR}n[l,l + 9] = 0}. 
On B{n) we can approximate Sn{t) by iS^e(t) + Z^^{t): 



sup \S^{t) - S^],{t) + Z^\t)\ < 6. 

iG[0,T] 



Further, let 7^*f^ = 

G2{n) 
G;{n) 



B = 

{S^'^.it) + (t) : t G M}. We define 



(5.39) 
(5.40) 



{7^5 n ((1 - 5, 1 + e) u (1 + e - 5, 1 + + 5)) ^ 0}, 
Gl«^{7^5^ (1,1 + ^^)^0}, 

Giinf n {7^*', n (1, 1 + ^^) = 0} = Gi(n)'= n G2(n)^ 



(5.41) 



These three events form a partition of the probability space. The reason for this partition 
is the following. If Gi{n) fl B{n) happens, then the process + Z^'^ intersects 5- 
neighbourhoods of 1 or 1 + and it is therefore not possible to decide if E{n) is true. 
On the other hand. 



E{ny 



(5.42) 



B{n)nG^{n) =^ E{n) and B{n)nG2{n) 

as can be seen from ()5.40|) . Therefore 

¥[B{n) n G3(r2)|T] < ¥[E{n)\T] < P[S(n)"|r] + F[Gi{n)\T] + F[Gs{n)\T]. (5.43) 

Thanks to Corollarv 15. 31 it is possible to estimate the probabilities of Gi{n) and G3{n). 
Indeed, + Z^'^ converge weakly in the Skorokhod topology to an a-stable subordi- 
nator. The probability that the subordinator hits any of the boundary points 1 ± 5, 
1 + 6* ± 5 is zero. Therefore, it follows from the weak convergence and Corollaries 15.31 
and IA.5I that 



lim P[G3(n)|T] = AsU 



1-5 

TeTs 



(5.44) 
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Similarly, by ()A.9|1 . 



P[Gi(n)|r] < 1 - AsU ([^) + 1 - AsU { [ll + 1 ) ^ ^^'^^^ 

Since 5 is arbitrary, the claim of the theorem follows from ()5.4H|1 - ()5.45|1 and from the 
continuity of AsIq,(-). This proves aging on the complete graph for k < 



5.2 The a-stable subordinator as a universal clock 

We will now give a general set of conditions which ensures that the result and the proof 
that we have given in the simple case of the complete graph apply. Consider an arbitrary 
sequence of graphs G„ and the Bouchaud trap models on them, BTM(G„, r„, 0). We set 
z/ = 1, therefore at every vertex x the Markov chain X waits an exponentially distributed 
time with mean r^/dr^ [dx is the degree of x) and then it jumps to one of the neighbouring 
vertices with an equal probability. We suppose that Assumption 12. Hl holds. in particular 
that means that are i.i.d. 

We want to prove aging for this model at a time scale = t^^n) using the same 
strategy as for the complete graph. That means, depending on the sequence Gn and the 
time scale 1^: 

• to divide traps into three groups: shallow, deep, and very deep, 

• to prove that the shallow traps can be ignored since the time spent there is negli- 
gible, 

• to prove that the very deep traps can be ignored because they are not visited in 
the proper time scale, 

• to show that the contribution of the deep traps to the time change can be ap- 
proximated by a stable subordinator, which will show the convergence of the clock 
process, 

• and finally, to deduce aging for the two-time function R from this convergence. 

The idea is that, as in the case of the complete graph, the time change Sn{j) should 
be dominated by a relatively small number of large contributions coming from the deep 
traps. 

On the other hand we want to stay as general as possible: we do not want to use 
any particular properties of the graph. Therefore, we formulate six conditions on the 
sequence G„ and the simple random walk on it. If these conditions are verified, 
the proof of aging can be finished in the spirit of Section IK. 1.41 Proving these conditions 
should be dependent on the graphs Gn- 

To formulate the conditions it is necessary to choose several objects that depend on 
the particular sequence Gn and on the observation time scale t^{n). 
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First, it is necessary to fix a (random) time C,n up to which we observe Yn- Second, a 
scale g{n) for deep traps should be chosen according to Gn and t^{n). This scale defines 
the set of the deep traps by 

Tf (n) := {xeVn-- eg{n) < r,. < Mg{n)}. (5.46) 

A possible generalisation of the definition ()5.46|) is described in the remark after the 
Theorem 15.71 

For the complete graph we used g{n) = n'^ and C,n ~ n'^°' (see the proof of Lemma IK3|l . 

As before, we use TMiji) = {x : r^. > Mg{n)} to denote the set of very deep traps. 
Similarly, we write T'^{n) = {x : < eg{n)} for the set of shallow traps. 

It should be possible to almost ignore these two sets. This is ensured by the following 
two conditions. Compare them with Lemmas 15.51 and 15.61 

Condition 1. There is a function h{e) satisfying lim£^o^(^) = 0, such that for a.e. re- 
alisation of T and for all n large enough 

r ^" 1 
E[5Ze,ry„(,)]l{y„(z) G T%n)} rj < /i(e)t^(n). (5.47) 

j=0 

In words, the expected time spent in the shallow traps before ^„ is small with respect 
to twin). 

Let if„(y4), A C V„, denotes the hitting time of A by the simple random walk y„, 

H^{A) ■= inf{i > : r„(z) G A}. (5.48) 

Condition 2. Given for any 5 > there exists M large enough such that for 
a.e. realisation of r and for all n large 

^[Hn{TM{n))<^T] <S. (5.49) 

We should now ensure that the contribution of the deep traps to the time change 
can be approximated by an a-stable subordinator. The following facts were crucial for 
the proof in the complete graph case: asymptotically, as n oo, 

• every time a deep trap is visited, its depth is independent of depths of the previ- 
ously visited deep traps. 

• the probability that a deep trap with depth ug{n) = un'^ is visited is proportional 
to 

• for fixed e and M only a finite number of deep traps was visited before the time- 
horizon ^. 
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However, these facts could not be true for a general graph. Here, after leaving a deep 
trap, the process Y typically hits this trap with larger probability than any other deep 
trap. On recurrent graphs, such as Z^, it even visits this trap a number of times that 
diverges with n. 

To overcome this problem we "group" the successive visits of one deep trap. The 
time spend during these visits there will then be considered as one contribution to the 
clock process. We define r„(j) as the sequence of times when a new deep trap is visited, 
r„(0) = 0, and 

r„(^) = min {j > r^{t - 1) : G Tf (n) \ {F„(r„(? - 1))}}. (5.50) 

We use (n to denote the largest j such that r„(j) < 

Cn ■■= max{j : r„(j) < (5.51) 

We define the process Un{j) that records the trajectory of F„ (and thus of X„) restricted 
to the deep traps, 

[/„(j) :=y;(r„(j)), jGNo. (5.52) 

Finally, define the score Sn{j) be the time that X„ spends at site Un{j) between steps 
r„(j) and r„(j + 1), 

rn{j+l) 

Snii) ■= ^ e,ry„(i)]l{F„(i) = f/„(j)}, j < Cn- (5.53) 

i=rn(j) 

Denoting by G^{x,y) the Green's function of the simple random walk Yn killed on the 
first visit to a set A C V^, it is easy to observe that s„(j) has the exponential distribution 
with mean 

^;'^t/n(i)G^M\{f/nO-)}(f^nO'), UnU))- (5.54) 

Since Conditionsnand|21ensure that the visits of deep traps determine the behaviour 
of the time change 5'„(j), the sum Yl\=i ^n{i) can be considered as a good approximation 
of Snivnij))- The next condition guarantees that the scores Sn{i) have a good asymptotic 
behaviour, i.e. are independent and have the right tail. 

Condition 3. Let {soo{i) : i G N) be an i.i.d. sequence given by Soo{i) '■= eia^{i), 
where cr^^{i) is a sequence of i.i.d. random variables taking values between e and M 
with common distribution function 

-a _ -a p-a _ -a 

and Cj is an i.i.d. sequence of exponential, mean-one random variables independent of 
. Then there exists a constant /C > such that for all £, M and for a.e. r, the sequence 
{sn{j)/'tw{n), j G N) converges as n — > cxd in law to the sequence (/Csoo(j), j G N). (For 
notational convenience we define s„(j) = Soo{j) for all j > C„.) 
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The last three conditions will ensure that the approximation by the a-stable subor- 
dinator is relevant for the aging. 



First, we need that 5'„(r„(C„)) is larger than {l + 9)tw{n). Since r„(C„) > YliZi 
and Sn{i) are easier to control than 5'„(r„(j)) we require 

Condition 4. For a.e. r and for any fixed > 0, 5 > it is possible to choose ^„ such 
that for all e small and M large enough, and for defined in ()5.5H) 



P| 

i=l 



[ J2 ^" W ^ + ^)'^M\r] >l-5. (5.56) 



Second, to prove aging for the two-point function R we need to show that for any 
time t' between 5'„(r„(j)) and Sn{rn{j + 1)) the probability that X„(t') = Un{j) is large. 
For a formal statement of this claim we need some definitions. Let t'^ be a deterministic 
time sequence satisfying tw{n)/2 < t'^ < (1 + 9)tw{n), and let 5 > 0. We define jn G N 
by 

SnirnUn)) < < Snij^Un + 1)) - St^. (5-57) 

and in = oo if ()5.57p is not satisfied for any integer. Let An{5) be an event defined by 

An{5) = {0 < 3n < Cn}. (5.58) 

We require 

Condition 5. For any 5 it is possible to choose e small and M large enough such that 
for a.e. r and all n large enough 

P[X„(4) = Un{jn)\An{5),T] > I - 5. (5.59) 

The last condition that we need excludes repetitions in the sequence Un- 
Condition 6. For any fixed e and M and a.e. r 

lim P[30 < i, J < Cn such that i ^ j and f/„(z) = ?7„(j)|t] = 0. (5.60) 



We have formulated the six conditions that are inspired by the complete graph proof. 
It should be then not surprising that they imply the same result as for the complete 
graph: 

Theorem 5.7 (Aging on general graphs). Assume that ConditionsU\^ holds. Then 
for a. e. realisation of the random environment r 

lim Rnit^n), (1 + 9)t^{ny, r) = AsU(l/l + 9). (5.61) 
This theorem can be prove d in a very similar way as Theorem 15.11 The complete 



proof can be found in |BC06a 
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Remark. 1. The set T^\n) as defined in ()5.4(jj) is sometimes too large and not aU 
conditions that we formulated can be verified easily. Typically this happens when two 
points of the top are too close to each other with non-negligible probability. In such 
case it is useful to define a set B{n) C T^'\n) of bad traps which cause the difficulties. 
If it is then possible to verify Conditions for the set T'^\n) := T^^{n) \ B{n) and, 
moreover, if B{n) satisfies for a.e. r a similar condition as the set of the very deep traps 

limsupP[//„(5(n)) < ^„|t] = 0, (5.62) 



then the conclusions of Theorem 15. 71 hold without change. 

2. The last pair of conditions is, in principal, necessary only for a "post-processing". 
If they are not verified, it is possible to prove aging for a top- dependent correlation 
function 

K{tn,,t^ + t-.T) = ^[^3 ■.Sn{r^U))<t^<t^ + t<Sn{r,^{3 + l))\T\, (5.63) 

which gives the probability that at most one site in the top is visited by X„ during the 
observed time interval. A two-point function similar to R!^ was considered in jBBG03aj . 



Theorem 5.8. If only ConditionsU^^hold, then T-a.s. 

lim K{t^{n), (1 + e)t^{ny, r) = AsU(l/l + 6). (5.64) 

n— +00 

5.3 Potential-theoretic characterisation 

In the previous paragraph we have stated six conditions that allow to prove aging on an 
arbitrary sequence of graphs. It is however not clear if these conditions can be verified for 
any concrete model. In fact, they are satisfied for the graphs: Gn = G = 7/ with c? > 1; 
for Gn a large torus in d dimensions, G„ = 7//nU^] for the n-dimensional hypercube, 
Gn = { — 1, 1}" among others (included of course the complete graph!). These examples 
will be developed in the next section, but we wont be able to give proofs, which can be 



pretty difficult (see jBCMn6( BC06aJ). Rather, we want to give here a few hints with 
potential-theoretic fiavor on how to verify our conditions. 

We want mainly discuss the crucial Conditional the convergence of the scores s„(z) 
to the i.i.d. sequence Soo- In the discussion we suppose that Gn are finite and sufficiently 
regular. Observe first that the set of deep traps T*^ is a random cloud on Vn (i-e. set of 
points chosen independently from V„, we assume that Tx are i.i.d.). The intensity p„ of 
this cloud depends on n; under Assumption 12.31 p„ ~ p^^g{n)^°' = (e^" — M~°')g{n)~°'. 
We, obviously, need that p„ 0, so that the random cloud is sufficiently sparse, but also 
Pn|V„| — >■ oo, so that the mean size of the random cloud diverges (otherwise fiuctuations 
of the random depths are important and no a.s. convergence holds). 

Going back to sequence Sn{i), we have already remarked that conditionally on Un{j), 
the j^^ visited deep trap, the score Sn{j) is exponentially distributed with mean 

''n^^ulijfUr.U)GTM(n)\{U^{j)}iUn{j), Un{j)), (5.65) 
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where G^{-, ■) is the Green's function of the simple random walk y„ that is killed on the 
first hit of the set A C V„. 

If the graphs Gn are "sufficiently regular" , then all vertices have similar degree, that 
means, e.g., that there is a scale d{n) such that dx/d{n) is uniformly bounded from 
and oo. Then by setting z/„ = d{n)~^ we can, at least theoretically, ignore the first two 
terms in This problem does not appear in all examples we consider, there always 

dx = const (n). Hence, it remains to control T[7^(j) and the Green's function in ()5.65|) . 

For the Green's function, one typically prove that there is a scale f{n) ~ tw{n) / g{n) 
such that T-a.s. 

/(n)-^G^M(„)\l^}(x,a;) const. (5.66) 

uniformly for all x £ T^{n). This result is again reasonable if the graph Gn is sufficiently 
regular and finite, and the cloud T^{n) is very diluted. 

To control the distribution of Tu^^(^j) consider first an arbitrary random cloud C V„ 
with intensity cp„. This random cloud will represent the set of the deep traps or its 
subsets. Recall that Hn{A) denotes the hitting time of A by Yn- Let P^^ be the law of 
Yn started at x. Suppose that it is possible to show for all m > and some scale r{n) 
independent of c 



sup 



X 



P. 



Hn{An \ {x}) ^ ■ 

r~\ — - ^ 

r[n) 



— exp(— cm) "^°°> 0, T — a.s., (5.67) 



so that the distribution of the normalised hitting time converges uniformly to the expo- 
nential distribution with mean c~^. This is again a reasonable property for very diluted 
clouds. If fl5.67p holds for all c then the lack-of-memory property of the exponential dis- 
tribution allows to prove the following claim: Let A^ be a random cloud with intensity 
(a + h)pn and let Bn C An be its sub-cloud with intensity hpn- Then uniformly for all 

X G An 

¥.[H{Bn \ {x}) < H{An \ {x})] (5.68) 
This claim yields, e.g. under Assumption 12.31 that for u G {e,M) 

^uAj-i)[ru.U)<u]^§^, (5.69) 

which is exactly what gives the "right tail" and the asymptotic independence of Sn{i)- 
The Condition IHl is then consequence of the claims of the last three paragraphs. 

If the graphs Gn are infinite one typically cannot prove uniformity in ()5.66|1 and 
()5.67|) . One can however restrict only to those deep traps that are reachable in ^„ steps 
from the starting position. On this set the uniform control is usually possible. 

There is also a heuristic reason for ConditionlHl As the random cloud T^{n) becomes 
more diluted, the hitting measure of this cloud charges more and more points because the 
random walk can "pass more easily around" . If this happens, the law-of-large-numbers- 
type arguments hold not only for the whole sets T^{n) but also for these sets as sampled 
by the hitting measure. That is why Condition El holds. 



39 



Remark. This is not true on Z! The simple random walk on Z cannot "pass around": 
it hits necessarily one of two points of the random cloud that are "neighbours" of its 
starting position. This explains the special properties of the BTM on Z and also the 
necessity of averaging in Theorems 13 . 1 1 1 and 13 . 1 31 

The Conditions 121 and E] are also easy consequences of claim ()5.67|) . One need first fix 
= mr{n) with m large enough to accumulate large enough (but finite) number K of 

"independent" scores Sn{i) in order to have with a large probability '^f^i Sn{i) > I + 0. 

This satisfies Condition 0] Then one fix M large enough, such that the intensity c{M)pn 

of TM{n) satisfies c(M) > 6^^m. Then Condition |21 holds. 
The preceding discussion can be summarised as follows: 

Claim. // ()5.66p and fl5.67|) can be checked for T^''^{n) on Gn, then, under Assump- 
tion \2.^A Conditions can he verified. 

The Condition n does not follow directly from ()5.(j7j) since Poisson clouds with larger 
intensity than p„ should be considered. The "slicing strategy" as presented in the proof 
for the complete graph however usually works. The remaining Conditions El and El are 
not substantial and we do not discuss them here. 

Remark also that Condition [H together with ()5.(ifi|l and ()5.fi7|l allows to prove the 
approximation of the clock process by an a-stable subordinator. This approximation is 
not a consequence of Conditions 1-6 only. 

There are at least two met hods of proof for facts ()5.66|) and ()5.67p . The first is the 



coarse-graining procedure of [ BCM06 that is explained in Section I^Til The advantage of 
this procedure is that it should work on many different graphs. It is however relatively 
technical and many special cases should be treated apart. O ne can also use the formula 



given by Matthews |Mat88j . This method is used in |BC06a to prove aging in the REM 



and on the torus in 1?. It however applies only if Yn "completely forgets" the position 
of Un{i) before hitting f/„(2 + 1), i.e. that the hitting measure of the random cloud is 
essentially uniform. Therefore this method does not apply e.g. on Z*^. 



6 Applications of the arcsine law 

We now describe two examples where the approach of Section El can be used to prove 
aging. 

6.1 Aging in the REM 

The Random Energy Model is the simplest mean-field model for spin-glasses and its 
static behaviour is well understood. The studies of dynamics are much more sparse. 
The first proof of aging in the REM was given in jK K(;n3airR"R(;n3bj . based on renewal 
theory. The general approach of Section El gives another, shorter proof. This approach 
allows to prove aging on a broader range of time scales, but on the other hand do not 
include quite exactly the results of |BBG03al IBBGOSb] . We will compare both results 
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later. Before doing it, let us define the model and give some motivation why and in 
which ranges of times and temperatures aging occurs. 

The Random Energy model is a mean-field model of a spin-glass. It consists of n spins 
that can take values —1 or 1, that is configurations of the REM are elements of V„ = 
{ — 1,1}". The energies {Ex,x G Vn} of the configurations are i.i.d. random variables. 
The standard choice of the marginal distribution of is centred normal distribution 
with variance n. We will, however, deviate from the standard choice to simplify the 
computations and we will assume that —E^ are i.i.d. positive random variables with the 
common distribution given by 

F[-EJ^/^ >u] = e""'/^ u>0. (6.1) 

We then define 

r, = exp(-/3E,). (6.2) 

The distribution ()6.1|) has almost the same tail behaviour as the normal distribution. 
As we already know, it is the tail behaviour of (and thus of E^j.) that is responsible 
for aging. Therefore, the use of "faked normal distribution" ()6.1|) is not substantial for 
our discussion. Remark also that a similar trick, i.e. to take E^ to be minimum of and 
the normal variable, was used in jBBGOSb] . 

For the dynamics of the REM we require that only one spin can be flipped at a given 
moment. This corresponds to 

n 

£n = {(x, y) e V^ 5^ |x, - y,\ = 2}, (6.3) 

i=l 

where (xi,...,a;„) are the values of individual spins. We use G„ to denote the n- 
dimensional hypercube (V„, £n)- There are many choices for the dynamics of REM, that 
has the Gibbs measure r as a reversible measure. We will naturally consider the trap 
model dynamics BTM(G„, r, 0), which is one of the simplest choices. We fix i/„ = 1/n 
in definition ()2.2j) . so that Tx is the mean waiting time at x. We always suppose that 

r„(0) = X„(0) = l = (l,...,l). (6.4) 

6.1.1 Short time scales 

Theorem 6.1. Let the parameters a G (0, 1) and (5 > be such that 

3/4 < a^/5V21og2 < 1. (6.5) 

Define 

tu,{n) := exp(a/5^n). (6.6) 

Then, for a.e. t, 

lim Rnit^n), (1 + e)tM) = AsU(l/l + 6). (6.7) 
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Let us first explain the appearance of scale tw{n) together with one problem that is 
specially related to REM. We have seen in Section El that aging occurs only if Tj. are 
sufficiently heavy-tailed. This certainly fails to be true for the REM: an easy calculation 



gives F[tx > u] 



u 



■ \ogu/2l3'^n 



which decreases faster than any polynomial. It is therefore 



clear that, if the system is given enough time to explore a large part of the configuration 
space and thus to discover the absence of heavy tails, then no aging occurs, at least not 
in our picture. On the other hand, at shorter time scales the system does not feel the 
non-existence of heavy tails as can be seen from the following estimate. Let a > 0, then 



exp < — 



log^ u 

2/3% 



logM + aP'^n' 



— a logyj 



and therefore 



P 



> U 



U 



oo 



(6.9) 



In view of the fact that the simple random walk on the hypercube almost never back- 
tracks, it seems reasonable to presume that if the process had time to make only ap- 



proximately e 



steps, then it has no time to discover the absence of heavy tails and 



aging could be observed. The above theorem shows this presumption to be true. 

Let us remark that there is much stronger relation between "random exponentials" 
Tx and heavy-tailed random variables. Let G N) be an i.i.d. sequence of centred 

normal random variables with variance one. It was proved in |BBM05j that for some 
properly chosen Z{n) and iV(n) the normalised sum 



N{n) 



Z n 



('n^ ^ 



(6.10) 



i=l 



converges as — > oo in law to an a-stable distribution with a depending on j3 and N{n). 
We show that the same is true for the properly normalised clock process S{n), which is a 
properly normalised sum of correlated random variables, more precisely of i.i.d. random 
variables sampled by a random walk. 

In view of fj6.9|) it is easy to fix objects for which Conditions HHHl should be verified: 
we define 



(n,a) 



exp(Q;/5^n), 
mexp(a^/3^n/2), 

{xGV„:r, e(£,M)e"^'"}. 



(6.11) 
(6.12) 

(6.13) 



The Theorem 16.11 is then the consequence of the following proposition and Theo- 
rem EIZI 

Proposition 6.2. Let a and (3 he as in Theorem \f).l\ Then for any 9 it is possible to 
choose m large enough such that ConditionsU^^ hold for F-a.e. r. 
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We believe that the range of the vahdity ()(j.5|l of the Theorem lfj.il is not the broadest 
possible. The upper bound 1 is correct. If a^/5^/21og2 > 1, then ^ 2". That means 
that the state space V„ becomes too small and the process can feel its finiteness. On the 
other hand, the lower-bound 3/4 is purely technical and can probably be improved. 

Observe also that the condition ()fi.5|l can be rewritten as 



-^i5c\Pm <I3< «"'^c, (6.14) 



where Pc = a/2 log 2 is the critical temperature of the usual REM. This, in particular, 
means that aging can be observed in REM also above the critical temperature, (3 < f3c- 

The proof of this proposition in BCn6a follows the strategy outlined in SectionEland 
uses the results of Matthews fMat88j for the fine control of the simple random walk on 
the hypercube. In particular, ()5.67j] is a consequence of the following potential-theoretic 
result which might be of independent interest. 

Proposition 6.3. (i) Let for all n > 1 sets An C Vn be such that \An\ = Pn'^"' with 
"densities" pn satisfying lim^^oo Pn2'''" = p G (0, oo) for some 7 G (1/2, 1). Let further 
the sets An satisfy the minimal distance condition 

min{d{x,y) : x,y E An} > (c^(7) + e)n (6.15) 

for some small constant e > and for the unique solution uj{'~f) of 

a;logcu + (l-cj)log(l-cj) + log2 = (27-l)log2, cjG (0,1/2). (6.16) 

Then for all s > 



lim max 

n^oo x^An 



E, 



exp ( - ^Hn{An \ {x}))] - I = 0. (6.17) 



That means that the hitting time Hn{An \ {x})/2'^" is asymptotically exponentially dis- 
tributed with mean 1/p. 

(a) If An are random clouds with intensity pn such that lim^^oo Pn2'^" = P ^ (0, C)o) 
and 7 G (3/4, 1), then the assumptions of (i) are a.s. satisfied. 

This result also explains the appearance of the lower bound 3/4 in the range of 
the validity of Theorem 16.11 the set of deep traps satisfies the assumptions of Proposi- 
tion IHHii) only if log 2 > 3/4. 



6.1.2 Long time scales 

The result of |BB(7n3a| IBB(T03a] deals with the longest possible time scales where aging 
appears in REM. The continuous-time Markov process X is replaced by a discrete-time 
process X', which at every step has the possibility not to move. The number of tries 
before leaving x has geometrical distribution with mean t^. As n 00, this dynamics 
differs very little from the usual trap model dynamics. The random environment is given 
by 

Tx = exp{|3^/nma.x{Ex, 0)), (6.18) 
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where are i.i.d. centred normal random variables with variance one. 
To define the relevant time and depth scales we set 

Un{E) = /3,v^ + E//3,v^- log(47rnlog2)/2/5,v^. (6.19) 

We define the set of deep traps {the top) by 

Tn{E) = {xe Vn : E, > Un{E)}. (6.20) 

The function Un{E) is chosen in such way that, as n — > cxd and E is kept fixed the 
distribution of \Tn{E)\ converges to the Poisson distribution with mean that depends 
only on E. The mean diverges ii E —oo afterwards. The initial position of the 
discrete-time Markov chain X' is chosen to be uniformly distributed in Tn{E). 

A different correlation function considered in |BBG03a] : If x„(fc) denotes the last 
trap visited by X' before step k then 

U'^{k,k + l,E;T) = r[{X'{i) : z G {A; + 1, . . . , A; + /} H (T„(E) \ {x„(A;)}) = 0|r]. 

(6.21) 

It is essentially the same as the function R' (see (I5.63|) ) 
This is the main aging result of |RR(;n3bj 



Theorem 6.4. For any (3 > (3c = V2 log 2 and any e > 

U'^{cJ,{1 + 9)cJ,E;t) 



lim lim lim P 



Asl/3.//3(l/l + 



> 6 



(6.22) 



where Cn ~ Ql^V^-^niE) ^ 



In fact one can see that the dynamics of the REM when observed only on the top 
Tn{E) can be approximated very well when n ^ oo and ii^ — — oo by a BTM on the 
complete graph with M = |T„(^)| ~ e"^ vertices (see |BEI02, BBGOSal). 

Let us now compare the results of Theorems 16. II and 16.41 First, different correlation 
functions R and 11' are considered. This difference is not substantial, we believe that it is 
possible to eliminate the top dependence (i.e. to convert something like R' to something 
like R) of ()6.22|) by some post-processing in the direction of Conditional 

The a.s. convergence in Theorem l6.1l is stronger than the convergence in probability in 
Theorem 16.41 It is a consequence of the fact that much larger time scales are considered 
and the set of the deep traps Tn{E) is finite for fixed E. Therefore, we cannot use 
law-of-large-numbers-type arguments for the number of deep traps with depth in a fixed 
interval. We have seen this effects already in the case of the complete graph (see ()5.8p ). 

The main difference between the two theorems is in the considered top sizes and 
time scales. In Theorem 16.41 the size of the top is kept bounded as n ^ oo. This allows 
to apply "lumping techniques" to describe the properties of the projection of a simple 
random walk on the hypercube to the top, that is to prove that that an equivalent of 
the process f/„ (see ()5.52|) ) converges to the simple random walk on the complete graph 
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with the vertex set Tn{E). In Theorem (or more precisely, in Proposition l(j.2j) the 
size of the top T^\n) increases exponentially with n. This makes the application of the 
lumping not possible and techniques based on Matthews results, i.e. Proposition 16.31 
should be used. 

The time scale c„ ~ e^^""^^) ~ ^ppcu+pE/Pc Theorem lfj.4l corresponds to the case 
/ (^c = 1 and is much larger than the scale tu,{n) ~ e"^^" = e^^^""^. These scales 
approach if a(3/f3c tends to 1, which is the upper limit of the validity of Theorem 16.11 
It would be possible to improve this theorem by setting ty^{n) = e^^^"'f{n) with some 
f{n) — i> as n — i> oo sufficiently fast, but even then tw{n) -C c^. Exactly at a(3 / (3c = 1 
Theorem 16.11 does not hold. As we have already remarked, in this case it is necessary to 
use the double-limit procedure and the convergence in probability. 

6.1.3 Open questions and conjectures 

The Theorems 16.11 and 16.41 give rigorous proofs of aging in the REM. They are however 
only partly satisfactory. It would be nice to replace the RHT dynamics by a more 
physical dynamics, like e.g. Glauber, or, at least, to explore the a 7^ case. We believe 
that the long-time behaviour of the model should not change dramatically, however we 
do not know any proof of it. The problem is that the Markov chain X becomes a time 
change of a random walk in random environment on the hypercube. Moreover, the clock 
process and the random walk are dependent. 

Another natural direction of research is to extend the results for the RHT dynamics 
on the REM to other mean-field spin-glasses, like the SK model or the p-spin SK model. 
In these models the energies of the spin configurations are no longer independent. 
We strongly believe that the approach of Section El can be applied, at least for large 
p and for well chosen time scales. These scales should be short enough not to feel the 
extreme values of the -E^'s which rule the (model-dependent) statics, but long enough 
for the convergence to a stable subordinator to take place for the clock-process. The 
difficulty is to verify Conditions 1-6 if the E^^s are not i.i.d. The assumption that the 
Eri-^s are independent is used twice in the proof for the REM. First, we use it to verify 
Condition ^ that is to prove that the time spent in the shallow traps is small. We 
believe that this condition stays valid also for dependent spin-glass models. The second 
use of the independence is more substantial. It is used to describe the geometrical 
structure of the set of the deep traps. More exactly it is used to bound from below the 
minimal distance between deep traps and to show that the number of traps in T„"+^(n) is 
proportional to It is an open question if these properties remain valid for dependent 
spin-glasses. 

6.2 Aging on large tori 

Another graph where the approach of Section El can be used to prove aging is a torus 
m Z'^. For convenience we will consider only d = 2 here, although similar results are 
expected to hold for li > 3. 
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Let Gn = iVni^n) be the two-dimensional torus of size 2^ with nearest-neighbours 
edges, i.e. V„ = 1? jl^l? , and edge (x, y) is in En iff 

2 

mod2" = 1. (6.23) 

i=\ 

We use ?/) to denote the graph distance of x, ?/ G V„. Let further r = {r"}, x G V^, 
n G N, be a collection of positive i.i.d. random variables satisfying Assumption 12.31 We 
consider the Bouchaud trap model, BTM(G'„, r, 0) with v = 1/4. 

Theorem 6.5 (Aging on the torus). Let t^{n) = 22"/"-f^i-{7/«) 7 g (0,1/6). 
Then for F-a.e. realisation of the random environment r 

lim R{t^{n), (1 + e)t^{ny, r) = AsU(l/l + 6). (6.24) 

n— »oo 

The theorem follows from the following proposition whose proof can be found in 
BCn6aj . 

Proposition 6.6. For any 9 there exist m large enough such that Conditions 1-6 hold 
for 

tJn) = 22'^/°ni-^/", = m22"n^-^, 

(6 25) 

(n) = {x G V„ : G {e, M)22'^/"n-^/"}, 

The main motivation for Theorem 16.51 was to extend the range of aging scales on 
and mainly to really explore the extreme values of the random landscape. Namely, the 
BTM on the whole lattice does not find the deepest traps that are close to its starting 
position. In the first 2^" steps, it gets to the distance 2" and visits 0(2^"/ log(2^")) = 
0(2^"/n) sites. Therefore, the deepest visited trap has a depth of order 2^"/°/n^/", 
which is much smaller that the depth of the deepest trap in the disk with radius 2", that 
is 2^"/°. Eventually, the process visits also this deepest trap, however it will be too late. 
This trap will no longer be relevant for the time change since much deeper traps will 
have already be visited. The deepest trap is relevant only if the random walk stays in 
the neighbourhood of its starting point a long enough time. One way to force it to stay 
is to change to the torus. By changing the size of the torus relatively to the number 
of considered steps, i.e. by changing 7, different depth scales become relevant for aging. 

The range of possible values of 7 G (0,1/6) has, as in the REM natural 
bound and an artificial one. It is natural that 7 < cannot be considered, since if the 
simple walk makes more than 2^" (log 2")^"*"^ steps, e > 0, inside the torus of size 2", its 
occupation probabilities are very close to the uniform measure on the torus, that is the 
process is almost in equilibrium. The other bound, 7 = 1/6, comes from the techniques 
that we use. We do not believe it to be meaningful since we expect the theorem to hold 
for all 7 > 0. Actually, the result for 7 > 1 follows easily from the Theorem 14.51 for 
the whole lattice. In this case the size of the torus is much larger than So that, the 
process has no time to discover the difference between the torus and Z^. We also know 
that Theorem 16.51 holds also in the window [1/6, 1] since it can be proved by the same 
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methods as for Z^, [BCMQ6, . Nevertheless the complete proof in this window has never 
been written. 

The 7 = case corresponds to the longest possible time scales. We expect that a 
similar result as Theorem 16.41 is valid. 

The proof of Proposition 16.61 uses again Matthews' results: the following equivalent 
of Proposition 16.31 can be proved. 

Proposition 6.7. (i) Let An C V„ be such that \An\ = Pn'^^^ with the density pn 
satisfying lim^^oo 2'^"'n~'^pn = p for some 7 G (0, 1/6) and p G (0, 00). Let further An 
satisfy the minimal distance condition 



mm{d{x,y) : x,y & An} > 2"'n 
for some /t > 0. Then, forIC = (2 log 2)^-*^, 

s 



lim max 

n— >oo xeAn 



exp 



22n^l-7 



j-HiAn \ {x}) 



ICp 



ICp 



0. 



(6.26) 



(6.27) 



(ii) If An are, in addition, random clouds with the densities given above, then the 
minimal distance condition is a.s. satisfied for all n large. 



Appendix. Subordinators 

We use frequently the theory of increasing Levy processes in these notes. We summarise 
in this appendix the facts that are important for us. For a complete treatment of this 
theory the reader is referred to the beautiful book by Bertoin |Ber96j . 

Definition A.l. We say that \^ is a Levy process if for every s,t > 0, the increment 
V{t + s) — V{t) is independent of the process {Viu), < u < t) and has the same law 
as V{s). That means in particular V{0) = 0. 

We work only with the class of increasing Levy processes, so called subordinators. 
There is a classical one-to-one correspondence between subordinators and the set of pairs 
(d, /i), where d > and p is a. measure on (0, 00), satisfying 



(1 A x)p{dx) < 00. (A.l) 

The law of a subordinator is uniquely determined by the Laplace transform of V{t), 

E[e-^^W] = e-**("), (A.2) 

where the Laplace exponent 

POO 

$(A)=d+/ (1 -e-^^)/i(dx). (A.S) 



47 



The constant d corresponds to the deterministic constant drift. All processes appearing 
in these notes have no drift, therefore we suppose always d = 0. The measure /i is called 
the Levy measure of the subordinator V. 

There are two important families of subordinators. The first consists of the stable 
subordinators. A subordinator is stable with index a G (0, 1) if for some c > its 
Laplace exponent satisfies 



m = cA" = — r / (1 - e-^-)x"^-- dx. (A.4) 

r(l - a) Jo 

Here F is the usual Gamma-function. 

The second important family of subordinators are the compound Poisson processes. 
They correspond to finite Levy measures, /i((0, oo)) < oo. In this case V can be 
constructed from a Poisson point process J on (0, cxo) with constant intensity : = 
/i((0, oo)) and a family of i.i.d. random variables Sj with marginal Z~^fi as follows. Let 
J = {xi, i G N}, xi < X2 < . . . , and xq = 0. Then V is the process with V^(0) = that 
is constant on all intervals (xj, Xj+i), and at Xi it jumps by Sj, i.e. V{xi) — V{xi—) = Sj. 

We need to deduce convergence of subordinators from the convergence of Levy mea- 
sures. 

Lemma A. 2. Let Vn be subordinators with Levy measures fin- Suppose that the sequence 
fin converges weakly to some measure fx satisfying ()A.1|) . Then Vn converge to V weakly 
in the Skorokhod topology on D = D{[0,T), M) for all final instants T > 0. 

Sketch of the proof. To check the convergence on the space of cadlag path D endowed 
with Skorokhod topology, it is necessary check two facts: (a) the convergence of finite- 
dimensional distributions, and (b) tightness. To check (a) it is sufficient to look at 
distributions at one fixed time, since Vn have independent, stationary increments. From 
the weak convergence of /i„ it follows that for all A > 0, 

]E[g-AV„(i)j ^ g-i/(l-e-^-)/.„(da;) "^°o, g-t / (l-e-^-)Mdx) ^ ]E[e--^^(*)]^ (A.5) 

which implies the weak convergence of Vn{t). Since Vn are increasing, to check the 
tightness it is sufficient to check the tightness of Vn(T), which is equivalent to 

lim lim E[e-^^"(^)] = 1. (A.6) 

This is easy to verify using the weak convergence of fin and the validity of (jA.lj) for /i. □ 

Definition A. 3 (The generalised arcsine distributions). For any a G (0,1), the 
generalised arcsine distribution with parameter a is the distribution on [0, 1] with density 

^.-^(1-.)-. (A.7) 

71 

We use AsIq, to denote its distribution function, 

sin rvTT 

AsU(m):=/ ^^M""^(l-M)-"dM, mG[0,1]. (A.8) 
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Note that AsIq(z) = tt ^ sin (0:77)5(2;; a,l — a) where B{z] a, h) is the incomplete Beta 
function. It is easy to see that 

AsU(2) = l-0((l-2)^-"), asz^l. (A.9) 

The following fact is crucial for us. 

Proposition A. 4 (The arcsine law). Let V be an a-stable subordinator and letT{x) = 
inf{t : V{t) > x}. Then the random variable V{T{x)—)/x has the generalised arcsine 
distribution with parameter a. 

This proposition has an important corollary. 

Corollary A. 5. The probability that the a-stable subordinator V jumps over interval 
[a, b] (i.e. there is no t & W such that V{t) G [a, b]) is equal to 

F[V{T{b)-) <a]= AsU(a/6). (A.IO) 

Sketch of the proof of Proposition \A.4\ Consider the potential measure U of the subor- 
dinator V, 



U{A)=E / l{V{t) e A}dt . (A.ll) 
Its Laplace transform is given by 

poo |- roo ^ 1 

/ e"^"t/(dx) = E / e-^^W dt = — — . (A.12) 
Jo Uo ^ *(A) 

Define further fi{x) = fi{{x, 00)). Then 



00 



e-'^/i(x)dx = $(A)/A. (A.13) 

Fix a; > 0. For every < y < x < z, we can write 

F[V{T{x)-) G dy, V{T{x)) G d^] = U{dy)ij{dz - y). (A.14) 

(For a proof of this intuitively obvious claim see p. 76 of |Ber96j .) Define now At{x) = 
x~^V{T{tx)~) and consider its "double" Laplace transform 

/■oo 

A{q,X)= e-«*E[exp(-A^(x))]dt. (A.15) 
Jo 

This Laplace transform can be explicitly calculated. Indeed using ()A.12|1 and ()A.13j) we 
obtain 

/■oo ptx 

A{q,X)= e-^' e-^y/''ji{tx-y)U{dy)dt 



Jo 

00 ^00 



JO 



$(g/x) (^{q/x) 



X $((A + g)/x) q/x g$((g + A)/x)' 
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Using that V is a-stable, i.e. = cx°', we get 

Observing that 7r~^ sin(Q;7r) = (r(Q;)r(l — a))~^ yields the claim of the proposition. □ 
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